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ABSTRACT 


Static var compensators (SVC) are traditionally placed at the generator terminals to 
mitigate subsynchronous oscillations (SSO) in series compensated electric power 
transmission systems. An auxiliary control [14] of the SVC employing generator rotor 
frequency feedback is shown to be adequate for damping SSO. 

It has been reported [15] that an SVC located at the midpoint of a series 
compensated line can be utilized for dual purposes of damping SSO and stability 
enhancement. Controllers of these SVCs have been designed based on a combination of 
line current signal and computed value of generator internal frequency (CIF). Both these 
signals utilize local measurements, as the current trend in power systems practice is to 
use local signals for the sake of reliability. 

In this thesis an altogether new concept of SVC control employing a remote 
generator frequency signal transmitted over a telecom line, is presented for damping of 
SSO in series compensated transmission systems. 

The IEEE First SSR Benchmark System is suitably modified to include an SVC at 
the midpoint of transmission line. Auxiliary controllers are designed for various signals 
such as line current, computed internal frequency (CIF) and remote frequency signal. 

Their effectiveness is then investigated for damping all the torsional modes at all critical 
levels of series compensation. Eigenvalue analysis is utilized to examine the system 
stability in each case. 

It is concluded in this thesis that this remote rotor frequency signal together with 
line current signal can be successfully applied for damping all the torsional modes at all 
the critical levels of series compensation. The performance of remote rotor frequency 
signal is much better than the computed internal frequency signal. This remote frequency 
signal is observed to be efficient over a reasonable range of telecom delays. 

This concept of remotely transmitted signal is being investigated for the first time 
in the control of static var compensators to damp subsynchronous oscillations. 
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Chapter 1 

INTRODUCTION 


1.1 General 


The generation and consumption of electrical energy has been increasing at a 
tremendous rate throughout the world during the past many decades. The 
increasing need of transmitting greater and greater amounts of power has led to 
an increase in the transmission voltages at EHV and UHV levels. Series 
compensation is a powerful tool to improve the performance of EHV lines. It 
consists of capacitors connected in series with the line at suitable locations. 

1.2 Advantages of Series Compensation 

The use of series compensation leads to a number of advantages as regards 
power transfer capability, voltage regulation, etc. 

(a) Increase in Transmission Capacity 

The power transfer Pi over a line is given by 


R = 


Ex E 2 


sin 8 


X \ 


( 1 1 ) 


where E^ZS = sending end voltage 
E 2 Z0 = receiving end voltage 
X e = line reactance 

8 = phase angle between £, and E 2 

If a series capacitor having capacitance X c is inserted, the net series 
reactance becomes X e -X c and the power transfer P 2 is given by 



From eqns.(l.l) and (1.2), it is seen that for the same magnitude of E x , 
E 2 and 8 , P 2 is higher than P x . The increase in power transfer is given by 


P L= X, ^ 1 
P x X e -X Q \-k 


d-3) 


The factor k equals — - and is known as the degree of series 

X e 

compensation. Studies on transmission systems have revealed that in view of 
stable system operation, protection and economics, the practical degree of series 
compensation is chosen in the range of 40 % to 70 %. 

(b) Improvement of System Stability 

From eqns. (1.1) and (1.2), it is seen that for the same amount of power 
transfer and the same values of E x and E 2 , the phase angle 8 in the case of 
series compensated line is less than that for the uncompensated line. It is 
established from equal area criterion [17] that a lower angle 8 implies better 
system stability. Some other methods for improvement of system stability are also 
used. These include reducing the reactance of generators and transformers, use of 
bundled conductors and increasing the number of circuits in parallel. Series 
compensation generally offers the most economic solution with respect to both the 
steady state and the transient stability. 

(c) Load Division between Parallel Circuits 

The reduction of series reactance by series compensation makes it a useful 
tool to balance the loading in parallel circuits. When a system is to be 
strengthened by the addition of a new line or when one of the existing circuits is 
to be adjusted for parallel operation in order to achieve maximum power transfer 
and minimize losses, series compensation can be used to a great advantage. It has 
been reported [18], that in Sweden, the cost of the first series compensation in 
the 420 kV system was entirely recovered due to the decreased losses in the 220 
kV system operating in parallel with the 420 kV system. 

1.3 Problems Associated with Series Compensation 

The use of series compensation introduces a few problems too . Some of these 
are : 

(a) Subsynchronous Resonance 

The series capacitor introduces a sub-synchronous frequency (proportional to 
the square root of the compensation) in the system. In some cases this frequency 
may interact with weak steam-turbine generator shafts and give rise to high 



torsional stresses. In hydro-turbine generators the risk of subsynchronous resonance 
is small because its natural frequency is lower than the range of frequencies in 
which subsynchronous torque is induced due to practical levels of series 
compensation. 

(b) Ferroresonan.ee 

When an unloaded or a lightly loaded transformer is energized through a 
series compensated line, ferroresonance may occur. The frequency of oscillation is 
an integral multiple of the system frequency [18] This can be suppressed by 
using shunt resistors across the capacitors or by short circuiting the capacitor 
temporarily through an isolator or a bypass breaker. 

(b) Line Protection 

Series compensation can lead to mal-operation of the distance relays of the 
line protection if the degree of compensation and capacitor location is not proper. 
To ensure correct operation of the distance protection the series compensation is 
limited to about 40 % ( when the installation is located at the middle of the line) 
and 30 % per bank (when the installations are located at one-third and two-third 
distances along the line ) 

(d) High Recovery Voltage 

Series capacitors produce high recovery voltages across the circuit breaker 
contacts. 

Inspite of the above problems series compensation is very widely used for 
EHV systems in the 400-500 kV range in many countries of the world. 

1.4 Subsynchronous Resonance and its Control 

The subsynchronous resonance (SSR) phenomena is usually associated with 
synchronous machines connected to series compensated transmission networks . It 
has been defined by IEEE SSR task force [1] as follows: 

" Subsynchronous resonance is an electric power system condition where the 
electric network exchanges energy with the turbine-generator at one or more of 
the natural frequencies of the combined system below the synchronous frequency 
of the system.” 

The definition includes any system condition that provides the opportunity 
for an exchange of energy at a given subsynchronous frequency. This includes 
what might be considered “natural” modes of oscillation that are due to the 
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inherent system characteristics, as well as "forced” modes of oscillation that are 
driven by a particular device or control system. The first SSR problem was 
experienced on December 9, 1970 resulting in the failure of a turbine generator 
shaft at the Mohave plant in southern California [12]. It was not until a second 
shaft failure occurred on October 26, 1971 that the real cause of the failure was 
recognized as subsynchronous resonance. 

The SSR phenomena is known to occur in two different forms • 

(a) Self excited or steady state SSR resulting from the following two 
mechanisms 

( i ) induction generator action 

(ii ) torsional interaction 

(b) Transient SSR 

To understand these different aspects let us consider a simple radial system 
consisting of a synchronous generator connected to an infinity bus through a 
series compensated transmission line, as shown in Fig. 1.1. 

The series resonant frequency f e of the electrical network is given by 


fe=fo 



fo 




where 


X s = reactance of generator and transformer 

X e = line reactance 

X L - reactance of series capacitor 

t 

X, - total inductive reactance = X s + X e 

X 

k = degree of series compensation = — - 

X „ 


f 0 = nominal system frequency 

It is to be noted that f e < f 0 since k < 1 . 


(1.3) 


(a) Steady State SSR 

Electrical subsynchronous currents of frequency f e flowing in the armature 
induce subsynchronous torques and currents in the rotor circuit having frequency 
f r given as 


fr=fo~fe 


(1.4) 



These rotor currents result in subsynchronous armature voltage components 
which may sustain or enhance the subsynchronous armature currents to produce 
the 'self-excitation’ effect or steady state SSR. Self excitation can be divided into 
two categories, one involving electrical system dynamics alone and the other 
involving both the electncal system and mechanical system (turbine-generator) 
dynamics. 

Induction Generator Action 

The induction generator effect relates only to the electrical system 
dynamics. This results from the apparent negative resistance characteristic of 
generators at frequencies below the system synchronous frequency At the system 
resonant frequency f e , if this apparent negative resistance exceeds the sum of the 
armature and network resistances, the subsynchronous armature currents may be 
negatively damped. 

Torsional Interaction 

This form of self excitation involves both the electrical and mechanical 
system dynamics. Torsional interaction problems may occur when the electrical 
resonant frequency f e is near the complement of a torsional resonant frequency 
/„ of the turbine-generator (T-G) shaft system. In this context, the word 
‘complement’ implies the difference between synchronous frequency and the 
torsional natural frequency. Under these conditions, a small voltage induced in the 
generator armature by rotor oscillations can result in a large subsynchronous 
current. When the net circuit resistance is positive, this current will produce a 
rotor torque which is phased to sustain or enhance the rotor oscillations. If the 
component of subsynchronous torque in phase with rotor velocity deviation equals 
or exceeds the inherent damping torque of the rotating system, the coupled 
electromechanical system will experience growing oscillations. 

It is important to note that the induction generator effect and torsional 
interaction are not mutually exclusive but will co-exist. They are treated separately 
for ease in analysis . 

(b) Transient SSR 

Transient SSR generally refers to transient torques on segments of the T-G 
shaft resulting from subsynchronous oscillating currents in the network caused by 
faults or switching operations. 

This usually occurs when the complement of the electrical network resonant 
frequency gets closely aligned with one of the torsional natural frequencies. 
Although, these transient torques decay with time their magnitudes are large and 
have the potential to cause serious damage to the T-G unit. 



Infact, even if torques such as those produced by steady state 
subsynchronous currents may be limited in magnitude and result in stresses 
within the elastic limit, failure of a shaft can be caused due to cyclic fatigue 
stress, which is cumulative. 

The phenomena of subsynchronous oscillations (SSO) is not restncted to 
series compensated systems alone. SSO has been reported to result on account of 
adverse interaction between the turbine-generator torsional system and HVDC 
control [19-21], dc converter loads and power system stabilizers [21]. An SSR 
condition can be created with shunt capacitor compensated systems if it is 
attempted to operate a long distance transmission system at a power angle in 
excess of 180° . 

1.4.1 General Countermeasures to SSR Problems : 

A wide variety of methods have been proposed in the literature [5] for 
counteracting SSR. These may be broadly classified as follows: 

(a) Static filter : This can take the form of a blocking filter (parallel-resonance 
type) in series with the generator, or damping circuits in parallel with the 
series capacitors. 

(b) Dynamic filter : This is an active device placed in series with the generator . 
It picks up a signal derived from rotor motion and produces a voltage in 
phase opposition so as to compensate for or even exceed the subsynchronous 
voltage generated in the armature. 

(c) Dynamic Stabilizer . This consists of thyristor modulated shunt reactors 
connected to the generator terminal Control of subsynchronous oscillations is 
achieved by modulating the thyristor switch firing angles, using signals derived 
from the generator shaft speed. 

(d) Excitation system damper : Generator excitation control is modulated by 
using a signal derived from the shaft speed so as to provide increased 
damping of torsional oscillations. 

(e) Protective relays -The. SSR condition is detected by a relay and the 
affected units are tripped. The relay may take several forms : detection of 
excessive torsional motion by sensing rotor speed or detection of SSR 
condition by sensing the armature current. 

(f) NGH scheme : This consists of a linear resistor in series with back-to-back 
thyristors connected across the series capacitor. The presence of SSR 
components in the capacitor voltage is detected and reduced by the action of 
the resistor discharge circuit. 



1.4.2 Control of SSR by SVC 

Static var compensators (SVCs) are shunt connected static reactive power 
generators and/or absorbers whose outputs are adjusted to exchange capacitive or 
inductive current in a continuous manner. This is done to maintain or control 
specific parameters of the electrical power system, typically bus voltage. The term 
"static” is used to indicate that SVCs, unlike synchronous compensators, have no 
moving or rotating main components. Thus an SVC consists of static var 
generator (SVG) or absorber devices and a suitable control device. 

A static var system (SVS) is a combination of different static and 
mechanically-switched var compensators whose outputs are coordinated. 

Static var compensators have been in use since the early 1960s. The SVC 
began being applied for transmission system voltage control in the late 1970s. 
Since that time, many SVCs have been applied worldwide for voltage control and, 
in some cases stability enhancement. The dynamic range of SVC is typically from 
60 MVAr to 600 MVAr. 

An SVC installed at the midpoint of a long transmission line can be 
utilized for damping of torsional oscillations in addition to voltage control [4], 
Different auxiliary control signals such as line current, computed internal 
frequency (CIF) of the generator [15], remote frequency signal and a suitable 
combination of them can be utilized to improve the damping of torsional 
oscillations significantly . 

1.4.3 Analytical Methods to study SSR 

Torsional interaction effects involve energy interchange between the 
generator shaft system and the inductances / capacitances of the network Therefore, 
the analysis of SSR problems requires representation of both the electromechanical 
dynamics of the generating units and the electromagnetic dynamics of the 
transmission network. 

Several methods have been used for the study of SSR [1,5]. A brief 
description of these methods is given below 

1.4.3.1 Eigenvalue (modal) Analysis 

Eigenvalue analysis uses the standard linear, state-space form of system 
equations and provides an appropriate tool for evaluating system conditions for the 
study of SSR, particularly for induction generator and torsional interaction effects. 

1. 4.3.2 Frequency Scanning 

This technique computes the equivalent impedance as seen from the 
internal buses of generators looking into the network, for different values of 
frequency. It gives information about the natural frequencies of the system and 
the tendency towards self-excitation and SSR. This approach is particularly suited 
for preliminary analysis of SSR problems. 



1. 4.3.3 Frequency response analysis of full system 

The stability of the subsynchronous oscillations is analyzed by using the 
multidimensional Nyquist criterion. This approach can handle detailed models and 
large systems: 

1.4.3. 4 Approximate frequency-domain analysis 

This rriethod analyzes the stability of individual torsional modes. It is 
limited to approximate detection of SSR. 

1.4.3. 5 Time domain analysis 

The electromagnetic transient program (EMTP) may be used to compute 
the transient time response. It allows very detailed representation of equipment, 
including nonlinear effects It is particularly suited for analyzing transient shaft 
torque due to SSR. However, this method is limited to the analysis of small 
systems. 

In this thesis, eigenvalue analysis is employed for the study of SSR, as 
this technique offers the following advantages : 

• Uses the state-space equations, making it possible to utilize many other 
analytical tools which can be applied in the same equation form 

• Computes all the exact modes of system oscillation in a single 
computation. 

• Can be arranged to perform a convenient parameter variation to study 
parameter sensitivities. 

• Can be used to plot root loci of eigenvalue movement in response to 
many different types of changes. 

1.5 Objectives and Scope of the Thesis 

The purpose of this thesis is to study the effect of SVC auxiliary 
controllers in damping all the torsional modes in the IEEE First SSR Benchmark 
System for all the critical levels of series compensation. 

While all auxiliary controllers reported in the literature [15,23] are based on 
locally available signals, in this thesis for the first time a new auxiliary signal is 
studied which is transmitted on telecom channels from a remote location. This 
signal is the frequency of the synchronous generator, and is termed remote 
frequency (RF) signal. 

In brief, the different objectives of this thesis are as follows : 

• To develop different component models of IEEE First SSR Benchmark 
system [2] for the study of SSR. 
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• To study the effect of a mid point located static var compensator (SVC) 
with voltage control alone on damping of torsional oscillations. 

• To examine the effectiveness of various local auxiliary control signals 
such as line current, computed internal frequency (CIF), remote generator 
frequency (RF) signal and a combination of them with respect to their 
ability to stabilize all the torsional modes. 

• To study the effect of telecom delays associated with remote frequency 
signal on the damping of torsional modes, rotor mode and electrical 
mode. 

• To study the effectiveness of auxiliary signals when SVC is placed at the 
generator end . 

1.6 Organization of Chapters 

Chapter 2 describes the development of models for the synchronous generator 
and its torsional system and the network for the IEEE First Benchmark Model for 
SSR analysis. The system eigenvalues with network tuned to different torsional 
modes are also presented . 

Chapter 3 presents the modeling of network, SVC voltage controller and 
auxiliary controller for IEEE First SSR Benchmark Model when SVC is placed at 
the mid point of transmission line. 

Chapter 4 studies the effect of different auxiliary signals such as line 
current, CIF, remote frequency (RF) signal and their combinations with respect to 
their ability to stabilize all the torsional modes. The effectiveness of auxiliary 
signals is examined at all the four critical series compensation levels for IEEE 
First SSR Benchmark system, i.e. 35.6 %, 54.8 %, 73 % and 90.8 %, when mode 4, 
mode 3, mode 2 and mode 1 respectively are critically destabilized. 

Chapter 5 studies the effect of telecom delays associated with the 
transmission of RF signal on the damping of various critical modes Effectiveness 
of auxiliary signals is also studied when SVC is placed at the generator end. 

A brief review of the various contributions of this thesis and the scope for 
future work has been consolidated in the sixth chapter. 


10 



Chapter 2 

SYSTEM MODELING 


2.1 Introduction 

This chapter presents a detailed formulation of the synchronous generator 
model, the multi-mass turbine system model and the network model for IEEE 
First Benchmark Model for SSR analysis. The overall model of the system is 
obtained by appropriately combining the above component models. The system 
model thus obtained is nonlinear in nature . 

2.2 Study System 

The study system is the IEEE First Benchmark Model for SSR analysis, 
depicted in Fig. 2.1. An 892.4 MVA synchronous generator is connected to an 
infinity bus via a highly compensated 500 kV transmission line. The mechanical 
system consists of a four stages steam turbine, the generator and a rotating exciter. 

The generator was represented by a Type - 59 synchronous machine source 
component. Two damper windings were provided in the q axis while at the d axis 
one damper and a field winding were considered. The mechanical system was 
represented by a multi mass spring dash pot system, with six lumped masses 
coupled by shaft sections of known torsional elasticity. Mechanical damping was 
assumed to be zero, to represent the worst damping condition 

2.3 Generator Model 

The mathematical model of the synchronous generator used here is as 
proposed in [6]. The schematic of the synchronous generator is shown in Fig 2.2. 
This shows three identical armature windings a, b, c and four rotor windings f, h, 
g and k. T coil represents the rotor field winding. The fictitious ‘g’ coil 
represents the effect of eddy currents which circulate in the solid steel of rotor. 

The ‘h’ and ‘k’ coils represent the damper winding along d and q axis 
respectively. All windings except the field winding are short circuited as they are 
not connected to voltage sources. The additional dummy coil ‘ c ’ has been 
introduced to handle transient saliency. The dummy coil is considered to be 
closed through a very high resistance and is assumed to have no mutual damping 
with other rotor coils. 



Turbine 



< eo 


hp : high pressure turbine lpl : low pressure turbine 

ip : intermediate pressure turbine lp2 : low pressure turbine 




The following assumptions are made : 


(i) Fundamental frequency mmf distribution is considered in the air gap . 

(ii) Saturation is ignored . 

Assumption (ii) can be relaxed but is generally used to simplify the 
analysis. The effects of machine damping and prime mover dynamics are small 
and can also be neglected for simplicity. Such a model is considered to be 
sufficiently accurate for the purpose of studying torsional interaction. 

The generator model described here can be further subdivided into • 

(a) stator circuits 

(b) rotor circuits 

(c) mechanical system 

2.3.1 Stator Circuits 


In this model, the stator of synchronous generator is represented by a 
dependent current source Is in parallel with an inductance Ls as shown in 
Fig. 2. 3. The dependent current source replaces the time varying coupling between 
the rotor windings and stator windings. It may be noted that Is is a (3x1) 
vector and Ls is a (3 x 3) matrix. These are expressed as 

Is=[/„ I, /,]'=/„ c + I,s (2.1) 

where 

c' = 7(2/3) [cos/9 cos(<9-2;r/3) cos(/9 + 2;r/3)] 
s' =7(2/3) [sin/9 sin(<9-2;r/3) sin(<9 + 2;r/3)] 


Id , I q are components of dependent current source along the d and q 
axis , respectively. 9 is the rotor angle. Subscript t indicates transpose . 


L s = 


L 


1 1 1 
1 1 1 
1 1 1 


+ - 


24/ 


2 

1 



( 2 . 2 ) 


Such a representation of machine can handle both the symmetrical and 
unsymmetrical networks equally well. If the external network connected to 
machine terminals is symmetrical, as considered in this case, a, b, c components 





R a 



(a) a Axis representation 


Ra 



(b) (3 Axis representation 

Ra 

| WV « — T 
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o 

(c) o Axis representation 

FIG. 2.4 EQUIVALENT REPRESENTATION OF 
SYNCHRONOUS MACHINE ON a, (3, 


AXIS 



can be transformed to a, (5 , o components. The advantage of this transformation 
is that all the three a,/3 , o component models are uncoupled. 

Moreover, a network is identical with /? network and is the same as 
positive sequence network model. 

Equivalent source representation of the machine on a,/3 , o axes is 
derived in [7] and shown in Fig. 2.4.The equivalent circuit consists of three 
meshes which are not mutually coupled. R u denotes the armature resistance. The 
currents in the meshes correspond to a ,j3 , o components of the armature 
currents. The relationship between a,j3 , o components and the phase currents 
i a , 4 and / c is given by 
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(2.3) 


The dependent current sources in frame of reference are defined as 

I a = I d cos0 + I q sin 6 (2.4) 

Ip = -IjSind + I q cosd (2.5) 

2.3.2 Rotor Circuits 

The rotor flux linkages are defined by 


'T / = a^ t +a^ h +b l v f +b 2 i J 
V h =a^ f+ a^ h +b,i d 

%=a s 'V„+ a< ,'V k +b s i, (2.6) 

T, = 0,^+0, V«+V, 

%=*,%+!>, i. 


where v f is the field excitation voltage . 

constants a, - a, , 6, - />,, are defined in Appendix A. The d and q 
components of the machine terminal current are defined by i d ,i g as 




i a cos 9 + i b cos(<9 - + i c cos (9 + ~) 

i a sin 6 + i h sin(<9 - + i c sin(<9 + 

3 3 

It is noted that currents i d and i q are defined with respect to the machine 

reference frame. However, to have a common axis of representation with the ac 
network these currents are transformed to D-Q frame of reference which is 
rotating at synchronous speed co 0 The following transformation is employed 


cos <5 -sinJ 
sin <7 cos 8 


(2.7) 


where i D , i Q are the respective components of machine current along D 
and Q axis . 8 is the angle by which d axis leads the D axis . 

Substituting eqn (2 7) in eqn. (2.6) and linearizing the resulting equation 
gives the state equation of the rotor circuits as 

Xr~ A.rXr^ B r\Ur\ Jc B riUri (2 8 ) 


where 

= ^ A % ]' • !&,=[** A H' 

W„ : = [ A ‘» A 'J ' 

Matrices Ar Br\'Br an< ^ Bri are defined in Appendix B,Sec.B.l. 

The output equation of the rotor subsystem is developed by utilizing the 
relationship between I d , I q and the rotor flux linkages as given in [6]. 


Ij +C 2^h 

I'^cP'+cM+cft 

where constants c, - c 5 are defined in Appendix A. 


(2.9) 


I d , I are now transformed to D-Q axis components I D and I Q , 
respectively, using the transformation 


V 


cos 5 

sin 8 



-sin£ 

cos 8 



( 2 . 10 ) 


The output equations are finally derived as eqns. (B.6), (B.7) and given by 


ia 



( 2 . 11 ) 

( 2 . 12 ) 


y l ?! Cr\Xr + Dr\Ur\ 


y R2 C R2 X R + D R2 U Rl + D R2 U R2 
where y RX = t A J ‘> A ! q) '-y R2 = [ AI o a/J ' 

Matrices C«, » C R2 ' D R r Dri » D R2 and Z) ff4 are defined m Appendix B, 
Sec. B.l . 


2.4 Mechanical System Model 

A study of the torsional interaction between turbine - generator and series 
compensated transmission network necessitates a detailed modeling of the 
mechanical shaft system. The mechanical system can be described using either the 
‘ multiresonant mechanical model (all -modes model)' [2,12] or the 'modal 
mechanical model’ [2,12]. A comparison of these models is given in [12,14], The 
all- modes model is utilized for studies presented in this chapter as it includes the 
effect of zeroeth mode and accounts for the coupling between different modes . 

In the multiresonant model the turbogenerator mechanical system is 
represented as a linear multi-mass spring dash pot system. Each major rotating 
element is modeled as a lumped mass represented by its inertia while every shaft 
element is modeled as a mass-less rotational spring with its stiffness expressed by 
the spring constant. Viscous damping of each mass and shaft segment is 
represented by a dash-pot damping A general turbine-generator shaft system 
model is shown in Fig. 2.5. 

The equation of motion of i th mass of the n-mass rotating system is given 
by 


where , 


M,s ,+ A- „(<* 8 ,-i) + D..8 ,+ A.,*. (8 r O + Z- 1 ,(8,~ 8 ,_,) 

+ K u yS l -S l+l ) = T, m -T ai =T l 

M, = moment of inertia of ith rotor . 


(2.13) 


T mi = mechanical torque applied on ith rotor. 

T a = electromagnetic torque on ith rotor due to electromagnetic energy 
conversion. For turbine rotors this is taken as zero . 

8 , = angle of rotation of ith mass with reference to a fixed reference 
frame 

7] =net applied torque on ith rotor. 

A , = self-damping coefficient of ith mass. 





D, j = Dj , - mutual damping coefficient between ith and jth masses . 

K, tJ = Kj , = spring constant of the shaft connecting ith and jth masses . 

A total of n equations similar to eqn.(2.13) can be written for n masses and 
combined to result in the following equation of motion 

MS + D8 + KS = T (2 14) 


where , 

M = diagonal matrix of the inertia constants associated with each mass 
D = tri-diagonal matrix of damping constants. 

K = tri-diagonal matrix of spring constants. 

8 = vector representing the angular position of each mass. 

T = vector representing the net applied torque on the rotors. 

Eqn. (2.14) which describes the torsional motion of masses can be 
linearized around a quiescent operating point and written in state space form as 

X u ~ A u X v/ Bm U\i (2-15) 


y i/ Cs uX\t 

where xu=i AS a 4 > AA \1' 

A Sq ,A 8 = incremental changes in angular displacement and velocity of the 

generator mass . respectively . 

u „ = Ar = [A7| AT,.... AT, AT,., AT.]' 

AT,=AT„,-AT„ 

It is noted that when turbine - governor dynamics is neglected 

A ^, = o 

while the above state space model has been described for a general n-mass 
mechanical system, the state and output equations for a specific 6-mass turbine- 
generator shaft system as utilized for the studies reported in this chapter, are 
derived in Appendix B,Sec.B.2. This mechanical system, as shown in Fig.B.l, 
relates to a high pressure turbine, an intermediate pressure turbine, two low 
pressure turbines, the generator and a rotating exciter. It is described by the 
following state space equations : 


Xm A.mX,u + Bm\Um\ + BmiUmi 


(2.16) 



Cm X u 

AS 2 A S 3 AS a AS. A S 6 a S l A S 2 AS 3 AA 4 AS 5 AS 6 ]' 

] ’ Uui = \ Al D ^ Z o] ’ y K f = [A^5 5 ] 

Ai d , Ai q are the incremental D-Q axis components of the current entering 

generator terminals in the series compensated system . 

A I D . Al () are the incremental D-Q axis components of the dependent 
current source 

Au • Bui’ Bui and (J u are defined in Appendix B , Sec .B.2 . 

2.5 Network Model 

The network model includes the model of generator stator, transformer, 
transmission line and the infinite bus. The generator transformer is represented by 
its leakage reactance. The magnetizing current is neglected. The transmission line 
is modeled by series R, L and C. The infinite bus is represented as a constant 
voltage and constant frequency source. 

Fig . 2.6 depicts the a- axis representation of the study system shown in 
Fig. 2.1. I a is the a component of dependent current source as described in the 
stator circuit model of synchronous machine . The current entering the generator is 
represented by i a . The terminal voltage at infinite bus is indicated by v la . The 
leakage inductance of the transformer at generator end is denoted by L n ,whereas 
L, 2 represents the inductance of reactor connected to infinite bus . 

The equations for the symmetrical network expressed on a- axis are 
written as 


Xu = l AS i 


r di a D . r " dl a 


C A = / 

d, - 


where L A = L + L n + L- n _ + L d 


R A = R + R a 


The above equations can be rewritten in matrix form as 

Xa^SlXCSlia+S^a 


(2.17) 

( 2 . 18 ) 


(2.19) 
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where Xa = [i a vj' 

Matrices £ } , and ^ are defined in Appendix B , Sec. B.3. 

Similarly, the equations for the p - network are 

X/3 = SlXp + S + Sl V lj3 ( 2 - 20 ) 

where Xp = [ b V J' 

Since the variables Xa » Xp are sinusoidal quantities in steady state, 

they result in a time varying system. In order to reduce the overall system to a 
time invariant one, it is necessary to transform the variables by a D - Q 
transformation on synchronously rotating reference frame . The D -Q components 
of the variable x are related to a - P components by the following transformation 

Xa J cos 6 r l sin6> r /l Xo (221) 

_X P \ |_ sin 8 r I cos^ r /J[xyJ 

where Xd = [i d v /J'’ Xq = [ 1 q v q\ 

I is an identity matrix of proper dimension. 8 , is the angle by which 
D-axis leads a - axis . 

Equations (2 19) and (2.20) are transformed to D-Q frame of reference 
using eqn(2.21) and then linearized. It is noted that Av !D = Av k; = 0 

The state equation of the series compensated network is finally obtained as 


X\ = A\X\ + BmUx\ + BsiUsi ( 2 - 22 ) 

where x . v = \x, > X c ]' ' Mv, = Vln Iq] 

Us 2 = [a4 -W.V, 

Matrices ^ , $ V1 and ft S2 are defined in Appendix B, Sec. B.3. 

The current at generator terminals constitutes the output of network model. 
The corresponding network output equation is given by 


y N =Csx 


s 


(2.23) 


^ A 



where y = [a/ 0 A i Q ] ' 

Matrix (J s is defined in Appendix B. Sec B.3. 


2.6 Derivation of System Model 


The interconnection between the various subsystems can be mathematically 
described by the following relationship 


Ur = Fry, ( 2 - 24 > 

where 

Ur = system mput vector = [ Ur] Uri Um Uxn Wm m J' 


= system input vector = j; 



Matrix F: is defined in Appendix B, Sec B.4 . 

The state and output equations of all the constituent subsystems are combined to 
give 

X/ = A,X, + B,Ur (225) 

y r =CiX, + D,u, ( 2 - 26 ) 

where Xr = \x R Xu xJ' 

Matrices Ai ' Bi ' Ci and Dr ^ defined in Appendix B, Sec . B.4 . 

Substituting eqn. (2.26) in eqn. (2.24) 

m,4-F,D,]"F,C,x, <2 - 27) 

Substituting eqn.(2.27) in eqn (2.25) results in 


X,= A X, 

where 

A = At + Bii I ~FiDi] x F rC r 


(2.28) 



ROTOR CIRCUITS 



INTERCONNECTION OF VARIOUS SUBSYSTEMS IN OVERALL SYSTEM MODEL 






2.7 Eigenvalue Analysis 


Eigenvalue analysis has been performed for IEEE First SSR Benchmark 
System. The operating conditions are chosen to be, generator operating at 803.2 
MW. There is no SVC present. When the level of series compensation is such 
that the electrical frequency of network coincides with the natural frequency of 
oscillation of any one of the torsional masses, the real part of eigenvalue 
corresponding to that mode becomes highly positive 

Table 2 1 presents the system eigenvalues when network is tuned to 
different Torsional Modes. 

TABLE 2.1 Eigenvalues of IEEE Benchmark System without SVC. 

P g =1.0 p.u. , PF=0.9 Lag 


Mode 

Identification 

Xc ( 4)=0.178 

p.u. 

X C (3)=0.274 

p.u. 

Xc(2)-0‘365 

p.u. 

Xc(ip0.454 

p.u. 

Mode5 

0.0 

0.6 

0.0 

0.0 


±j298.1767 

±j298.1767 

±j298.1767 

±j298.1767 

Mode4 

+1.4279 

+0.0095 

-0.0002 

-0.0023 


±j202.8760 

±j202.7796 

±j202.8626 

±j202.8954 

Mode3 

+0.0147 

+ 1.3816 

+0.0149 

+0.0010 


±j 160.7544 

±j 160.6349 

±j 160.3997 

±j 160.4985 

Mode2 

+0.0012 

+0.0072 

+0 6715 

+0.0064 


±j 127 0498 

±j 127.0920 

±j 127.0063 

±j 126.9361 

Model 

+0 0021 

+0.0190 

+0.1248 

+4 8860 


±j99.4066 

±j 99.6600 

±j 100.38 18 

+J99.2286 

ModeO 

-0 3321 

-0.3905 

-0.4799 

-0.6222 


±j 10.41 74 

±j 11.2879 

±j 12.3397 

±jl3.7131 

Electrical Mode 

-5.3254 

-4 9425 

-3.8999 

-7.3889 


±j202.7614 

±j 160.7794 

±j 126.9643 

±j99 2066 

Supersyn . Mode 

-4.7601 

-4.7973 

-4.8233 

-4.8436 


±j55 1.0233 

±j592.9337 

±j626.2365 

±j654.9739 

Generator 

-40.7906 

-40.9539 

-41.1552 

-41.4114 

Rotor Circuits 

-25.4118 

-25.4264 

-25.4457 

-25.4733 


-3.0801 

-3.3957 

-3.7929 

-4.3401 


-0.1287 

-0.1333 

-0.1429 

-0.1599 


-0.0676 

-0.0871 

-0.1016 

-0.1100 


From the above results, we see that IEEE First SSR Benchmark System 
is characterizes by 5 torsional modes corresponding to 5 masses of the torsional 
system. These torsional modes are called as model, mode 2, mode 3, mode 4 and 




mode 5 whose frequencies (in radian/sec) are 99, 127, 160, 202 and 298, 
respectively. There is one mode called as electrical mode, whose frequency 
matches exactly with that of critically tuned mode. Electrical mode results from 
the interaction between line inductance and series capacitor. 

The eigenvalue results presented in Table 2.1 are in close agreement with 
the published results [14]. 

2.8 CONCLUSIONS 

In this chapter, system model is developed for IEEE First SSR Benchmark 
system. The potential instabilities of this system are assessed by an eigenvalue 
program. This analysis shows that the system is characterized by four unstable 
torsional modes distributed over a wide frequency range. 

The eigenvalue results of Table 2.1 agree with the published results. 
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Chapter 3 


MODELING OF 
SVC CONTROL SYSTEM 


3.1 Introduction 

It has been reported earlier [15,23] that a midpoint located SVC can not 
only improve system stability but also damp subsynchronous oscillations in a 
single machine infinite bus system. However, this has been demonstrated with 
auxiliary controllers of SVC utilizing locally available signals and a computed 
generator frequency signal which has been shown to be quite effective. 

In this chapter a new concept is investigated where instead of computing 
the generator frequency, the actual rotor frequency signal is measured and 
transmitted over a telecom channel to the remotely located SVC. 

This chapter describes the modeling of SVC compensated network, SVC 
voltage control system and modeling of various SVC auxiliary controllers. 

3.2 Study System 

The study system is the IEEE First Benchmark System with an SVC 
placed at the mid point of transmission line . as depicted in Fig 3.1. 

3.3 System Model 

The system model is derived from the models of constituent subsystems : 

(i) Synchronous generator 

(a) Stator circuits 

(b) Rotor circuits 

(c) Mechanical system 

(ii) Network 

(iii) SVC 

(a) Voltage controller 

(b) Auxiliary feedback controller 
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Component models of the stator and rotor circuits of synchronous generator 
and mechanical system including representation of various turbine masses and their 
dynamics have already been derived in Chapter 2. The development of SVC 
compensated network, SVC voltage control system including TCR transients and 
all firing control and measurement delays have been descnbed in section 3.3.1 
and 3.3.2 respectively. 

Modeling of SVC auxiliary controller utilizing various auxiliary signals has 
been described in section 3.3.4 . 

3.3.1 Network Model 

The ac network considered for the study system of Fig. 3.1, is similar to 
that described in Chapter 2 except that the transmission line is segregated into two 
halves and an SVC is placed at the midpoint of line. The equivalent a - network 
is shown in Fig 3.2. 

R a denotes the armature resistance of the synchronous generator. I a 
is the a component of dependent current source as given in the model of 
synchronous machine. The current entering the generator is represented by i u . The 
voltage and current at infinite bus are represented by v ia and i [a respectively. L n 

and L n represent the leakage inductance of the transformer at generator end and 
inductance of reactor connected to infinite bus respectively. The SVC bus voltage 
is indicated by v 3tr and its current by L a . The two series capacitor voltages are 
represented by v, a and v 4a . The transmission line current entering SVC bus from 
the generator end is denoted by i 4a . The symbols R 1 ,Z 1 ,C, d and R 2 .L 2 ,C se2 
represents the resistance , inductance and capacitance of the two sections of 
transmission line 


The following differential equations are derived for the a - network 



R-, . 





(3.1) 



R 4 . 1 

-i. + — v. 


L, 


L* 


3 a 



L a dt 


(3.2) 


dv 2a _ 1 . _J_. __1_- 
dt C n ' a C n ha C n ha 


(3.3) 



dt 




‘I a 


(3.4) 





1 


(3.5) 


±la _ 

dt 



i 


a 


where L 4 = L n + + L” d 

R A =R\+R a > 


Lg Z/j*2 Ln 


c=c 


K 


The above equations can be written in matrix form as 


X a = S\X a + Slha + S'Ja + Si V \ a (3 6) 

where Xo = L L A* v i a V J' 

Matrices Si , S2 , S3 and S 4 are defined in Appendix C , Sec. C.l. 
Similarly, the equations for (3 - network are given by 

Xp-SiXfi + Sihp + SAe + S* v iji ( 3J ) 

where Xp = [hp h V v v V> V J' 

It is to be noted that vectors x^ , xp and matrices Si, S2 and S3 as 
given above are different from the similar named vectors and matrices as defined 
in chapter 2 corresponding to the system without SVC. 

Eqns.(3.6) and (3.7) are transformed to the synchronously rotating D -Q 
frame of reference using the transformation described by eqn. (2.21) and are 
subsequently linearized. Since the infinite bus voltage is constant. Avid , Aviq = 0 

The state equation of the series compensated network is finally derived as 


X.\ = As Xs + Bsi U VI + Bsi U V 2 + B V3 W.V3 


(3.8) 


where X v = [ 

Ad 

m d 

^ V 3/) ^ V 2D ^ V 4 I) Q 

Us 1 ~ 

Abo 

A Hq 

1 ' 

XI V2 ~ 

Vo 

A/y 

/ 

Un 3 ~ 

A 

A Iq. 

“ Usi 



Matrices As'Bni'Bni and J$ Ni are defined in Appendix C, Sec. C.l . 


The output variables of network model are the current at generator 
terminals and the SVC bus voltage. The corresponding network output equations 
are derived as eqns. (C.1,C.2) and given by 



y v j C- m Xs 

y N 2 ~ C \ 2 X X 

where J/ v] = K AzJ' and J / v2 = [av 30 Av 3 J' 

Matrices (J X] and (J %2 are defined in Appendix C, Sec. C.l 

3.3.2 Static Var Compensator 

SVC is a dynamic reactive power device which provides rapid and 
continuously controllable lagging and leading MVARs . Variation of reactive power 
can be achieved through both the active and passive control . 

3.3.2. 1 Operating Characteristic of SVC 

Fig. 3.3 shows the steady state control characteristic of an SVC, which 
represents the relationship between SVC current (or reactive power output) and the 
bus voltage. The steady state operating point is established at the intersection of 
the control characteristic with the system load line. For small deviations in the 
bus voltage around the controller set point V ref , the SVC current is regulated 

within the control range to provide inductive compensation for voltage rises and 
capacitive compensation for voltage drops. Thus SVC maintains the terminal 
voltage at the preset reference V nf . Large changes in the bus voltage force the 
SVC beyond its control range For severe undervoltages the SVC gets tranformed 
into an equivalent fixed capacitor having susceptance B [ - B c , while for large 
overvoltages, it reduces to an equivalent shunt reactor of susceptance B LKLiX - B c . 

The slope of the control characteristic essentially represents a compromise 
between SVC rating and the voltage stabilizing requirement. It also improves the 
current sharing between the SVCs operating in parallel. The slope of the control 
characteristic is typically chosen in the range of 1-10%. 

3.3.2.2 Modeling of SVC Voltage Controller 

A small signal model of a general SVC control system is depicted in Fig. 
3.5. The terminal voltage perturbation AV and the SVC incremental current 
weighted by the factor K D representing current droop are fed to the reference 

junction. T M represents the measurement time constant which for simplicity is 

assumed to be equal for both voltage and current measurements. The voltage 
regulator is assumed to be a proportional-integral (PI) controller. Thyristor control 

action is represented by an average dead time T D and a firing delay time T s . The 

variation in TCR susceptance is represented by A B . 


(3.9) 

(3.10) 




P : operating point 

B c = SVC capacitive susceptance 

B, = SVC inductive susceptance 


FIG. 3.3 STEADY STATE CONTROL 
CHARACTERISTIC OF SVC 



In detailed modeling of SVC, TCR transients are being included . The delays 
associated with S V S controller and measurement unit is also incorporated . The 
a ,P axis currents entering TCR from the network are expressed as 


L v 


dha 

dt 


+ R S ha = V 3 « 


(3.11) 


U 


d /, 


dt 


+ R S hp ~ V 3/3 


where L s , R s represent the inductance and resistance of the TCR , respectively . 

Transforming equation (3.11) to D-Q frame of reference using equation 
(2.21) and linearizing, gives 
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where Q = Quality Factor of TCR =- < % — , B = 


R s cOqL^ 


The following equations can be written for SYS control system 


where 


=1 = A Ke/~ Z 2 

A V rcf is the reference voltage perturbation 


(3.13) 


Av 3 - K d A/ 3 = z 2 + T u z 2 (3-14) 

where Av 3 , A ,/ 3 are the incremental magnitudes of SVS bus voltage and TCR 
current. 
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Substituting eqn. (3.13) in eqn. (3.16) gives 
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FIG . 3.4 GENERAL CONTROL SYSTEM BLOCK DIAGRAM FOR THYRISTORIZED SVC 















z 3 = AB + T d AB 


(3.18) 


or 


AB = — z 3 - — AB 
T T 

1 D 1 D 


(3.19) 


The magnitude of TCR current is given by 


h ~ V (bn + ho ) 


(3.20) 


Linearizing eqn. (3.20) 


a/ 3 =^ a, 


*30 


bo 


30 


(3.21) 


The SVS bus voltage v 3 is also expressed in terms of its D , Q axis 
components as 


V 3 =V( V 3/> 2+V 34? 2 ) 

Linearizing equation (3.22) gives 


(3.22) 


Av3= Im Av 3fi +^_ Av 3 £ 


(3.23) 


K 30 


•'30 


Eqns. (3.21) , (3.23) are substituted in eqn (3.15) to give the state equation 
corresponding to z 2 . 

The state and output equations of the SVC model can then be written as 

Xs = AsXs + BsxUv + Bsi U S 2 + Bs2 U S3 ( 3 ‘ 24 ) 

(3.25) 

r y V A * / V.V JLS A 

r 1 

where , 


+ DsUs\ 



Xs = [ A/ 3/J 

A h Q 2, % 2 Z 3 

ab]' 

Al bJ' 



A ^/ » W .V3 

= Av f , y s = [a/ 30 

O) 

. rn 

<1 


Matrices , g sl , J$ n , Qs and Ds are defined in Appendix C , Sec.C.2 



3.3.2.3 Choice of SVC Rating 

The dynamic range of SVC is chosen on the basis of reactive power 
requirement at the SVC bus to control voltage under steady state conditions. This 
information is obtained from load flow studies which, in addition, provide the 
voltage magnitudes and phase angles at various buses needed to compute the 
initial conditions in the system. The SVC bus is assumed to be a PV bus with 
zero power injection (SVC losses are neglected). 


Load flow is then conducted and SVC bus reactive power is computed for 
varying generator power outputs over whole range of series compensation from 
0 % to 90.8 % (a series compensation level corresponding to critical excitation of 
mode 1). The results are presented in Table 3.1. It is seen that as the generator 
real power increases from 100 MW to 803.2 MW the SVC reactive power output 
varies from -4.5 MVA inductive to 250.5 MVA capacitive. Hence, the dynamic 
range of SVC is chosen as 25 MVA inductive to 275 MVA capacitive. It is 
important to note that as the SVC configuration is FC-TCR, the reactor is rated 
larger than the fixed capacitance in order to provide net lagging VARs 

Table 3.1 

Requirement of Reactive Power at SVC Bus (in MVARs) 

K=% series compensation 


Level of Power 

Transfer 

(in MW) 

K=0.0 

K=35.6 

K=54.8 

K=73.0 

K=90.8 

100 

3.3 

1.7 

0.6 

-1.1 

-4.5 

500 

94.9 

67.0 

50.5 

31.9 

5.3 

803.2 

250.5 

179.3 

139.2 

96.8 

42.2 


The slope of the steady state V-I characteristic of SVC which is defined on 
nominal voltage base and the larger of the two reactive ratings, is taken as 1 % 
at 275 MVA capacitive reactive power output of SVC. 

3.4 Derivation of System Model 

The interconnection between the various subsystems can be mathematically 
described by the following relationship 


Ur F \y T 


(3.26) 





FIG .3.6 INTERCONNECTION OF VARIOUS SUBSYSTEMS IN OVERALL SYSTEM MODEL 







where 


si Ur2 Mm i Mm 2 Us'i Un2 Uni 


U 7 = system input vector = 
y f = system input vector = [y y„ y u y m y n y) 
Matrix p is defined in Appendix C ,Sec..C.3 . 


Us\ 


The state and output equations of ail the constituent subsystems are combined to 
give 


x? Ai x T + B iU t + B us2 + Busi 


y i CiXt+DiUt 


where 



Xm 


Xn 



(3.27) 

(3.28) 


Matrices At'Bt'B Ct 311(1 Dr are defined in Appendix C,Sec. C.3 . 
Substituting eqn. (3.28) in eqn . (3.26) 

U,=['-FtD,Y'FtCtXt <3-29) 

Substituting eqn.(3.29) in eqn. (3.27) results in 

X 7 = A X 7 + BUs2 + BUsi ( 3 - 3 °) 

where 

A = A,+B,[>-FrDr\''FrC r 

U S2 = Av W ’ U S2 = Av r 

where Av F is the output of the auxiliary controller 


3.5 SVC Auxiliary Controller Model 

This section deals with the development of a state space model for the 
SVC auxiliary controller. Fig. 3.7 depicts the block diagram of an SVC control system 
involving feedback of a general auxiliary signal u c through a controller transfer 
function G(s). In this analysis the perturbation in SVC reference voltage V nf is 
assumed to be zero. The state and output equations are derived for this controller 
with different auxiliary signals. Each signal is mathematically expressed as a 
function of the state variables of other subsystems. The various signals considered 



arc line current and CIF, which are locally available and generator frequency 
available at remote located SVC. 

3.5.1 Line Current Auxiliary Controller 

The magnitude of transmission line current z 4 entering SVC bus from the 
generator end, as shown in Fig. 3.1 is given by 

*4 =^w+<Mi (3.31) 

where z 4/) ,i 4y are the components of line current z 4 along D.Q axis, 
respectively . 

The auxiliary control signal chosen in this case is the perturbation in line 
current magnitude A z 4 , which is obtained by linearizing eqn.(3.31) 

M 4 =iML Aiw+ !m.Ai (3.32) 

*40 *40 

this signal can be expressed as 

u c = F ( rXr + FcmX a/ + FcnXn + FcsXs (333) 

where u c = A / 4 , fat = Fa, = Fes = 0 

Matrix J? CN is defined in Appendix C, Sec. C.4. 1 . 

The auxiliary controller is, a priori, assumed to be a simple. first order 
transfer function for this auxiliary signal and all the others described subsequently 
in this section. 


G(s) = 


A Vj. 

4 f s 


= K b ( 


1+57; 

1 + sT^ 


■) 


This controller can be equivalently expressed as 


(3.34) 


*2 


i zUlL. ) 

1 + sT 2 


(3.35) 


A block diagram of this auxiliary controller is displayed in Fig. 3.8. The 
corresponding state and output equations are given by 












1 + sT] 
1 + sT 2 








Xc = Ac Xc + Bc u c 


(3.36) 


y ( =CcX ( +Dc u ( ( 3 . 37 ) 

where Xc = M • >V = Av V 

Matrices Ac’Bc’Cc ^ Dc ^ defined in Appendix C, Sec.C.4.1 

3.5.2 Computed Internal Frequency (CIF) Auxiliary Controller 

Computed internal frequency is a control signal derived from the generating 
source frequency. Since the generating station and SVC are located far apart from 
each other, the CIF signal is computed from parameters which are available at the 
SVC bus. The quantities utilized for the computation of CIF signal are bus 
voltage, transmission line current at the SVC bus and the reactance between 
generator internal voltage and SVC terminals. The derivation procedure is given 
below. 


The a axis network of the study system including synchronous generator, 
transformers and transmission network shown in Fig. 3.2 is simplified to a form 
depicted in Fig. 3.9. For computing CIF signal, only the section between generator 
and SVC is considered. The resistances of the generator stator and transmission 
line is neglected. The fixed capacitor is considered as placed on the right of TCR 
and hence, is ignored in this analysis. The dependent current source representing 
the generator is transformed to an equivalent voltage source behind subtransient 
inductance. L E represents the total inductance between SVC and the equivalent 
voltage source and is given by 

L k =L” + L n +L (3.38) 

where 

L n d = subtransient inductance of the generator. 

Lj | = inductance of transformer at generator end. 

L— inductance of transmission line segment between the generator 
transformer and SVC . 

The or,/? axis components of the internal voltage e ] of synchronous generator is 
given by 


\a 

l/> 


= V +L ^ 
/: dt 

r di 4 p 

= ’ v L ' ~JT 


(3.39) 

(3.40) 



where v 3 „ , v 3fi are the components of SVC bus voltage v 3 along a ,J3 axis 

respectively. i Aa . i Ap are the components of line current / 4 along a ,/? axis 
respectively 


Transforming eqns. (3.39) and (3.40) to D-Q frame of reference using eqn. 
(2.21) results in 

e w = v 3 /> + L E i AD + co 0 L E i A Q (3.41) 

e \Q = v :-g + L h i A Q — O) 0 L r i Ai:) (3.42) 


where e U) , e UJ are the D and Q axis components of the internal voltage 


e, respectively. 


The angle S , of the internal voltage e, is computed as 

J,= tan-'(-^-) 

e \n 


(3.43) 


Linearizing eqn. (3.43) 


AS, 


Lmo. Ag --^-Ae 

2 ^ 1 Q 2 D 

e io e io 


(3.44) 


The computed internal frequency (CIF) signal A co l is then obtained by 
differentiating eqn. (3.44) 

e 

q ~ 


Ary,= — AS ,= — 2 
1 dt 1 e 2 


— Ae, 0 — Ae, 


w 


(3.45) 


10 


-10 


The expression for A ft), in general matrix notation is derived in Appendix 
C as eqn. (C.ll) and given by 


Uc = FcrXr + FcmXk{ + FcnX» + FcsXs 
where u c = A co , 

Matrices Fcr ’ Fat and Fat defmed in Appendix C, Sec. C..4.4 . 

The state and output eqn. Of the corresponding auxiliary controller are 
given by eqns.(3.36) and (3.37) respectively . 

3.5.3 Remote Frequency (RF) Auxiliary Controller 

The remote frequency signal is the actual measured value of generator 
rotor frequency which is transmitted over a telecom channel to the midpoint 
located SVC. This signal is associated with a time delay T td which includes both 
the measurement and telecom delay. Fig 3.10. represents the block diagram for 
remote frequency auxiliary controller . 
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The following eqns . can be written for remote frequency auxiliary 
controller . 


)^i +(^r— )Aty 5 
*TD 


z s ~ Ur ] z 5 + [B c ]^i 

Av f = [C r ]z 5 + [Z) r ]c| 


(3.47) 

(3-48) 

(3.49) 


The state and output equations are derived for this auxiliary controller as : 


Xc = A’cX c + B’c»c ( 3 - 50 ) 

yc = C'cXc ( 3 - 51 ) 

where X( = l z s e \\‘ » Tc = Av r > u c =Ato s 


Ac 


t _ 


Ac Be 

0 

1 1D J 


> Bc = 


0 

(^-) 


TD 


c\-\Cc al 

Matrices Ac , Be , Cc and Dc have been defined in Appendix C , Sec . C 4.1 


3.5.4 Composite Line Current and RF Auxiliary Controller 

The auxiliary signal in this case comprises a combination of line current 
and remote frequency (RF) signals with the objective of utilizing the beneficial 
contribution of both towards the damping of different torsional modes. 

The control scheme for composite controller is illustrated in block diagram 
form in Fig 3.11, where u n and u C2 correspond, respectively to the line current 
and RF control signals. Although, the auxiliary controller configurations are same 
for both the signals, they are referred as G(s) and G '(s) for line current and RF 
signals, respectively. The state and output equations of the two individual 
controllers, as derived in section 3.5.1, can be rewritten as 


Line Current 


Xn = AcXn + Bc u n 


o . * * 

I ^ 

(3.52) 
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>'n =CcXa + Dc u « 


(3.53) 


where % ( , — [z 51 ] , u n — A/ 4 , y n -Av F , 


Matrices Ac'Bc'Cc ^ Dc are defined in Sec. 3.3.1. 
Remote Frequency (RF) 

X<2 = Ac Xc2 + Be U < 2 

>’(2 =C( Xd + D( U C2 


(3.54) 

(3.55) 


where Xd = ^ z A ’ u c 2 ~^ 0) s » 3’c2 — ^ v f 2 

The state equation of the composite auxiliary controller is obtained by 
combining eqns.(3.52) and (3.54) 


t t 

x ( =Ac x ( + Bc Uc 

where x ( = lxn XcJ ’ Uc = i u n “ol 


(3.56) 
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Ac. 



Be 0 
o Be 


The auxiliary controller output is given by 


Av^ = Av n + Av F2 

The output equation of composite auxiliary controller is then obtained by 
adding eqns . (3.53) and (3.55) 


f * 

y c - =Cc Xc + De Uc 


(3.58) 


where y c - Av F 

C,’=tc< C<1 • Dc=\Dc Dc\ 
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I he modeling of composite line current and CIF signals is also derived in 
a similar manner. 

3.6 Derivation of System Model (With SVC Auxiliary 
Controller) 

The procedure given in Sec. 3.4 is utilized to combine the state and output 
equations of the series compensated network, SVC voltage controller and SVC 
auxiliary controller as developed in this chapter with those of the generator rotor 
circuit and mechanical system derived in Chapter 2, to result in the state equation 
of the overall system as 


x = Ax (3.59) 

where x = \x R Xu X N X s X r J' 

The general auxiliary controller eqn.(3.36) is simply appended to the 
system matrix A defined by eqn.(3.30) , to result in the modified system matrix. 

3.7 Conclusions 

In this chapter, the model for IEEE First Benchmark system together with 
an SVC placed at the midpoint of transmission line is developed. Further the 
models of SVC voltage controller and different auxiliary controllers viz. line 
current, computed internal frequency (CIF) and remote frequency (RF) signal are 
obtained. Models are also derived for composite auxiliary signals in the 
combination of line current and CIF, and line current and remote frequency. 
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Chapter 4 


SYSTEM STUDIES 

4.1 Introduction 


This chapter presents eigenvalue results for IEEE First SSR Benchmark 
System with an SVC placed at the midpoint of transmission line. SVC is 
employed in the Benchmark system with the objective to damp torsional 
oscillations It is important to note that if an SVC is placed at the midpoint of 
transmission line, then it also serves the purpose of improving dynamic stability. 
Eigenvalue results are presented for all the four critical series compensation levels 
for Benchmark system, i.e. 35.6 %, 54.8 %, 73 % and 90.8 %, at which network is 
tuned to mode 4, mode 3, mode 2 and mode 1, respectively. System studies are 
done with SVC voltage controller alone and with different auxiliary stabilizing 
signals viz . , line current, CIF, remote frequency signal and combination of line 
current with CIF and remote frequency signal, individually. 

4.2 Case Study 

The study system is the IEEE First SSR Benchmark System with an SVC 
placed at the midpoint of transmission line. The midpoint connected FC-TCR 
static var compensator is rated at 275 MVA leading to 25 MVA lagging . The 
operating value of generator power is taken to be 803 MW. The load flow data is 
given in Appendix D. 

The modeling of synchronous generator and turbine generator mechanical 
system has already been described in Chapter 2, whereas the modeling of network, 
SVC voltage controller and auxiliary controllers is described in chapter 3. The data 
corresponding to generator, rotor spring mass parameters and SVC voltage 
controller including all control system delays is given in Appendix D. To 
represent the worst damping condition all self and mutual damping coefficients 
are assumed to be zero . 

The objective of this case study is to examine the effectiveness of different 
auxiliary signals i.e. line current, CIF, remote frequency signal and combination of 
line current with CIF and remote frequency signal in damping torsional oscillation 
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modes. Eigenvalue analysis is employed to obtain the loci of critical system 
modes with varying controller parameters each taken one at a time, for each 
auxiliary controller. Based on these root loci, ranges are determined for the 
controller parameters ( Kb , Tj and T 2 ) which result in damping of all or most of 
the critical system modes. A suitable set of parameters lying in this range is then 
selected. The controller parameters which result in maximum damping of all or 
most of the critical system modes are chosen as optimal for the particular 
auxiliary controller. 


System studies are done at all the four critical series compensation levels : 
35.6%, 54.8%, 73% and 90.8% where torsional modes 4, 3, 2 and 1, 
respectively, are critically destabilized . 

4.2.1 Selection of voltage controller parameters 

For selecting voltage controller parameters Ki (integral gain) and K p 
(proportional gain ) the following technique is adopted. For selecting Ki , first the 
torsional system is removed from the rest of the system. The proportional gain K p 
is set to zero and Ki is varied over a wide range. From eigenvalue results root 
loci is plotted for sensitive eigenvalues. The value of Ki for which the rotor 
mode is most stable and all other system modes are also reasonably stable is 
chosen as optimal value of Ki For selecting K p , torsional system is incorporated 
and Ki is fixed at this optimal value. Then K p is varied and corresponding root 
loci is plotted. The value of K p at which the undamping of torsional modes is 
reduced is chosen as the optimal value of K p . These values ofK p and Ki are 
utilized for all the system studies . 

For IEEE First SSR Benchmark System with an SVC placed at the 
midpoint of transmission line, the optimal values of K p and K[ are found to be 
+2.6 and +90.0 respectively. 

4.2.2 Case 1 : Line Compensation = 35.6 % 

(Network tuned to mode 4 ) 

To begin with, the effect of SVC with pure voltage control in damping the 
torsional modes, is examined. Table 4.1 lists the system eigenvalues without SVC 
and with an SVC having voltage controller alone. It is seen that voltage controller 
alone has very little effect on the stability of torsional oscillations. 

Different auxiliary signals are now investigated in respect of their ability to 
stabilize the various torsional modes, these are : 

(i) line current 

(ii) CIF 

(iii) Remote frequency (RF) 

(iv) line current + CIF 

(v) line current + RF 
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4.2.2.1 Line Current Auxiliary Controller 


Eigenvalue analysis is performed to examine system stability as affected by 
the line current auxiliary controller parameters. Loci of eigenvalues sensitive to 
variation in K B , Ti and T 2 , each taken at a time, are plotted in Figs 4.1, 4.2 and 
4.3, respectively. 

From Fig. 4.1 , it is observed that for -0.41 <K B <-0.39 , all the torsional 
modes except mode 4 are successfully damped . As K B is reduced from -0.41 , 
the mode 4 tends towards stability, but the electrical mode tends towards 
instability. At K B = -0.42 rotor mode becomes unstable whereas at K B = -0.46 mode 
4 becomes stable, but this results in a high degree of destabilization of electrical 
mode. To study the effect of parameters T| and T 2 , a value of K B =-0.41 is 
selected at which all other modes except mode 4 are stable. 

The loci of critical eigenvalues asTi is varied is illustrated in Fig. 4 2. It is 
observed that for values of 0.096 < Ti <0.099, all other modes except mode 4 are 
stable. As T| is increased from 0.099, mode 4 tends towards stability but 
electrical mode and rotor mode tends towards instability. At Ti=0.12, mode 4 
becomes stable but rotor mode and electrical mode becomes unstable . 

Fig . 4.3 shows the behavior of critical roots with variation in T 2 . For 
0.0207 < T 2 <0.022, all other modes except mode 4 are stable. As T 2 is increased 
from 0.015 to 0.022, rotor mode and electrical mode tend towards stability, 
whereas mode 4 towards instability. The eigenvalue results corresponding to line 
current controller parameters which result in reduced undampings of mode 4 are 
shown in Table 4.1. 

From the above root loci, it is evident that line current controller is unable 
to damp all the system modes simultaneously at 35.6 % compensation. 

4.2.2.2 CIF Auxiliary Controller 

The variation of critical eigenvalues with variation in K B for CIF auxiliary 
controller is exhibited in Fig. 4.4. It is observed that as K B is increased from 
0.0025 to 0.0085, rotor mode and SVC mode tend towards instability whereas 
mode 2, mode 3, mode 4 and electrical mode tends towards stability . 

Fig. 4.5 shows the behavior of critical roots with variation in Tj.lt is seen 
that as T| is increased from 0.0001 to 0.0015, the rotor mode, and SVC mode 
tend towards instability whereas model, mode 2, mode 3, mode 4 and electrical 
mode tend towards stability . 

Fig . 4.6 shows the behavior of critical roots with variation in T 2 . It is seen 
that CIF could not stabilize modes 3, 4, SVC mode and rotor mode.System 
eigenvalues for a set of parameters chosen from the ranges described above are 
listed in Table 4.1 . 

4.2.23 Remote Frequency (RF) Auxiliary Controller 

The variation of critical eigenvalues with K B for RF auxiliary controller is 
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FIG . 4.7 ROOT LOCI WITH VARIATION IN K B FOR RF AUXILIARY 

CONTROLLER . LINE COMPENSATION = 35.6 % . T, = 0.043 Sec . , T 2 =9.0 E - 6 Sec 



59 






shown in Fig. 4.7. It is observed that for 0.007 <K B <0.011, all the torsional 
modes except mode 4 are successfully damped. As K B is reduced from K B = 0.007, 
modes 2 and 3 become unstable. 

Fig. 4.8 shows the behavior of critical roots with variation in T[. It is seen 
that all the modes except mode 4 are stable for 0.039<Ti <0.045. 

Fig. 4.9 shows the behavior of critical roots with variation in T 2 . It is seen 
that all the modes except mode 4 are stable for 9.0E-7<T 2 <9.0 E -5. 

System eigenvalues for a set of parameters chosen from the ranges described 
above are listed in table 4.1. 

For this signal, the telecom delay has been chosen to be 5.0 ms, and the 
same value has been utilized for all studies reported in this chapter. 


4.2.2.4 Line Current + CIF Auxiliary Controller 

A study of the eigenvalues given in Table 4.1 reveals that none of the 
auxiliary controller line current, CIF and RF alone can stabilize mode 4, if 
stability of all other modes are also to be ensured simultaneously. Therefore , the 
effect of composite line current and CIF signals is studied . 


The variation of critical eigenvalues with K' B for (LC +CIF ) auxiliary 
controller is shown in Fig. 4.10. It is observed that for 0.0080<K B <0.0085, all 
the torsional modes including mode 4 are successfully damped. As K' B is reduced 
from K'b = 0.0080, mode 4 becomes unstable. As K' B is increased from K' B - 
0.0085, the rotor mode becomes unstable. 

Fig. 4.11 shows the behavior of critical roots with variation in T'i. It is 

seen that all the system modes are damped for 0.0006<T' t <0.0007. AsT i is 
reduced from T'i =0.0006, electrical mode becomes unstable, whereas if T i is 
increased from 0.0007 electrical mode becomes unstable. 

Similarly, from the root loci displayed in Fig. 4.12, it is seen that whole 
system is stable for 0.018<T' 2 <0.020. As T' 2 is reduced from T' 2 -0.018, rotor 
mode becomes unstable, whereas if T' 2 is increased from 0.020 electrical mode 

becomes unstable. , „ . 

System eigenvalues for a set of parameters chosen from the ranges 

described above which ensure complete system stability are listed in 
Table 4.1. 

4.2.2.5 Line Current + Remote Frequ e ncy (RF) Auxiliary Contro ller 

The variation of critical eigenvalues with K ' B for ( ^ C ^ } the 

controller is shown in Fig. 4.13. It is observed that or . - B ^ red J ced 

torsional modes including mode 4 are successfu y ampe . B , 

from K' b = 0.014, mode 4 becomes unstable, however as K'b is increased from K B 

= 0.018, electrical mode becomes unstable. 
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i 4.9 ROOT LOCI WITH VARIATION IN T 2 FOR RF AUXILIARY 

CONTROLLER. LINE COMPENSATION = 35.6 %. K B = 0.009 , T, = 0.043 Sec 
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IG .4.10 ROOT LOCI WITH VARIATION IN K B FOR COMPOSITE 
(LINE CURRENT + CIF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION =35.6 % .T t ' = 0.0007 Sec . , T 2 =0.018 Sec. 
K n = -0.475 , T,- =0.084 Sec . , T 2 =0.0245 Sec . 
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FIG. 4. 11 ROOT LOCI WITH VARIATION IN T,' FOR COMPOSITE 

(LINE CURRENT + CIF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION = 35.6 %. K B ' = 0.0085 , T,' = 0.018 Sec. 
K B = -0.475 , T,- =0.084 Sec . , T 2 =0.0245 Sec . 



FIG. 4. 12 ROOT LOCI WITH VARIATION IN T 2 FOR COMPOSITE 

(LINE CURRENT + CIF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION = 35.6 % . K B =0.0085 , T,' = 0.0007 Sec . 
K b = -0.475 , Ti =0.084 Sec . , T 2 =0.0245 Sec . 
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ROOT LOCI WITH VARIATION IN K B FOR COMPOSITE 
(LINE CURRENT + RF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION = 35.6%.T,' = 0.009 Sec ., T 2 ' = 0.007 Sec ., T D =0.005 Sec 
K,j = -0.44 , T, =0.08 Sec. , T 2 =0.017 Sec. 
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■IG . 4.14 


ROOT LOCI WITH VARIATION IN Ti FOR COMPOSITE 
(LINE CURRENT + RF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION = 35.6 % . K B = 0.016 , T 2 - 0.007 Sec . , T D 0.005 Sec 
K b = -0.44 , Tr =0.08 Sec. , T 2 =0.017 Sec. 








•'IG .4.15 ROOT LOCI WITH VARIATION IN T 2 ' FOR COMPOSITE 
(LINE CURRENT + RF ) AUXILIARY CONTROLLER . 

LINE COMPENSATION = 35.6%.K B ' = 0.016 , T,’ = 0.009 Sec ., T D =0.005 Sec . 
K B = -0.44 , T r =0.08 Sec. , T 2 =0.017 Sec. 



FIG. 4. 16 ROOT LOCI WITH VARIATION IN K B FOR LINE CURRENT AUXILIARY 

CONTROLLER. LINE COMPENSATION = 90.8 %. Ti = 0.019 Sec T 2 = 0.017 Sec 
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SYSTEM EIGENVALUES FOR SVC VOLTAGE CONTROLLER AND DIFFERENT AUXILIARY 
CONTROLLERS. LINE SERIES COMPENSATION = 35.6 % 
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SYSTEM EIGENVALUES FOR SVC VOLTAGE CONTROLLER AND DIFFERENT AUXILIARY 

CONTROLLERS . LINE SERIES COMPENSATION = 54.8 % 
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Elect . Mode : Electrical Mode , S . S. Mode : Super-synchronous Mode 
Time Constants T| and T 2 are expressed in seconds 
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Parameters 

Elect . Mode : Electrical Mode , S S. Mode : Super-synchronous Mode 
Time Constants Ti and T 2 are expressed in seconds . 







Fig. 4.14 shows the behavior of critical roots with variation in T', It is 
seen that all the system modes are damped only for T', =0.009. AsT', is reduced 

from T', =0.009, mode 4 becomes unstable, whereas if T', is increased from 0.009 
electrical mode becomes unstable. 

Similarly, from the root loci displayed in Fig. 4.15, it is seen that whole 
system is stable only for V 2 =0.007. As T' 2 is reduced from T' 2 =0.007, electrical 
mode becomes unstable, whereas if T' 2 is increased from 0.007 mode 4 becomes 
unstable. System eigenvalues for a set of parameters chosen from the ranges 
described above are listed in Table 4.1. 

4.2.3 Case 2 : Line Compensation = 54.8 % 

(Network tuned to mode 3 ) 

Root loci studies were conducted for all signals as described in the 
previous section, and the best controllers were obtained. 

Eigenvalue results with these different auxiliary controllers are presented in 
Table 4.2. 

4.2.4 Case 3 : Line Compensation = 73.0 % 

(Network tuned to mode 2 ) 

Eigenvalue results with different optimized auxiliary controllers are 
presented in Table 4.3. 

4.2.5 Case 4 : Line Compensation = 90.8 % 

(Network tuned to mode 1 ) 

The variation of critical eigenvalues with K B for line current auxiliary 
controller is shown in Fig. 4.16. It is observed that for -1.20<K B <-1.05, all the 
system modes except electrical mode are successfully damped. As K B is increased 
from K b = -1.05, electrical mode becomes stable but mode 1 and mode 2 becomes 
unstable. 

Fig. 4.17 shows the behavior of critical roots with variation in Tj. It is 
seen that all the modes except mode 1 are stable for 0.01 9 <T[ <0.021 whereas for 
0.01 6 <T| <0.018 all other modes except electrical mode are damped. . 

Fig. 4.18 shows the behavior of critical roots with variation in T 2 . It is 
seen that for 0.0172<T 2 <.0.0174 all the system modes are stable. 

Thus , we see that at critical compensation level of 90.8 % , line current 
auxiliary controller alone is able to damp all the system modes . 

Eigenvalue results with line current auxiliary controller are presented in 

Table 4.4 . 

To sum up, the system eigenvalues for the stabilizing composite ( LC + RF ) 
auxiliary controller for different critical levels of series compensation are presented 
in Table 4.5. 


68 



k 

400 

Pts 

• T, 


1 

0.016 

300 

2 

0.017 


3 

0.018 

200 

4 

0.019 


5 

0.020 

175 

6 

0.021 

150 

125 


3 1 

100 


— O O 

6 



-50 -10 -1 - 0.1 


- 0.01 - 0.001 0 0.001 0.01 
Red) 


0.1 1 10 50 


FIG. 4. 17 ROOT LOCI WITH VARIATION IN T, FOR LINE CURRENT AUXILIARY 
CONTROLLER. LINE COMPENSATION =90.8 %. K B = -1.05 , T 2 = 0.017 Sec 


7 a*o 1 


-50 -10 -l 




7 a**' 


13 4 5 


L 

400 

Pts . 

t 2 


1 

0.0163 

300 

7 

0.0168 


8 

0.0170 

200 

9 

0.0172 

10 

0.0174 

175 

11 

0.0176 


12 

0.0180 

150 



125 

6 

7 

100 3 

l 


— v 

10 



0 JL- 

1 L 

1 1 1 !* 

A 1 1 !A CA 


FIG. 4. 18 ROOT LOCI WITH VARIATION IN ^2 p OR LINE CURRENT AUXILIARY 

CONTROLLER. LINE COMPENSATION =90.8 % . K B --1.05 , T, 0.019 Sec 





SYSTEM EIGENVALUES FOR SVC VOLTAGE CONTROLLER AND LINE CURRENT AUXILIARY 

CONTROLLER .LINE SERIES COMPENSATION = 90.8 % 
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TABLE 4.5 

SYSTEM EIGENVALUES FOR THE STABILIZING COMPOSITE (LC+RF) AUXILIARY CONTROLLERS FOR 

DIFFERENT CRITICAL LEVELS OF SERIES COMPENSATION 
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4.3 Discussions 


Mode 4 stabilization 


• From eigenvalue results presented in Table 4.1, for series compensation 
level of 35.6 % when mode 4 is critically destabilized, we see that SVC 
voltage controller alone has very little effect on the stability of torsional 
modes. This necessitates the need of auxiliary controller in addition to 
voltage controller. Therefore different auxiliary signals such as line 
current (LC), CIF, remote frequency (RF) and combination of line 
current with CIF and remote frequency are considered. It is observed 
that line current alone is unable to damp mode 4 , and CIF alone is 
unable to damp modes 3, 4, rotor mode and SVC mode while 
maintaining stability of all other modes. Remote frequency signal alone 
can not stabilize the electrical mode but enhances the stability of 
mode 4. Since none of the auxiliary signals alone are able to stabilize 
all the system modes, a need is felt to study the composite effect of 
two auxiliary signals. Hence the influence of (LC + CIF) and (LC +RF ) 
composite auxiliary controllers is examined. It is found that both (LC + 
CIF) and (LC + RF ) auxiliary controllers are able to stabilize the 
system. However, with (LC + CIF ) auxiliary controller the stability of 
mode 4 and rotor mode are lower in comparison to those ensured by 
(LC + RF ) auxiliary controller. 

Hence (LC + RF ) auxiliary controller is better than (LC + CIF ) auxiliary 
controller . 

Stabilization of mode 3 and mode 2 

• From eigenvalue results presented in Tables 4.2 and 4.3, for critical 
series compensation levels of 54.8 % and 73 % respectively, it is 
observed that line current alone is unable to stabilize the critically 
tuned mode while stabilizing the rest of the system modes.. It is seen 
that (LC + CIF ) combination is unable to stabilize the critically tuned 
mode at 54.8 % and 73 %, even though it could stabilize the system at 
35.6%. However, the (LC + RF) composite signal is able to stabilize 
all the system modes at 54.8 % and 73 % series compensation. 

Stabilization of mode 4 


• From eigenvalue results presented in Table 4.4, for critical series 

compensation level of 90.8 % when mode 1 is critically destabilized, it is 
seen that line current auxiliary controller alone is able to stabilize the 
system. 

As torsional oscillations are invariably caused by the interaction between 
the transmission line current and the turbine - generator mechanical 
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system, it can be expected that the line current signal will contain 
information regarding the torsional oscillations and hence constitutes a 
useful control signal . 

• From the eigen value results presented in Table 4.5, we see that over 
whole range of series compensation level, line current signal needs lead 
compensation (Ti > T 2 ) with negative gain constant (Kb). With increase 
in series compensation level, magnitude of gain increases and the phase 
lead (ratio of T 1 /T 2 ) decreases. 

• CIF signal needs lag compensation (Ti > T 2 ) for all cases, however gain 
constant is positive for mode 4 but negative with modes 3 and 2. 

Remote frequency (RF) signal needs lead compensation with positive 
gain constant. With increase in series compensation level, magnitude of 
gain decreases, whereas ratio of (T 1 /T 2 ) increases. 

4.4 CONCLUSIONS 

Eigenvalue analysis is conducted for the IEEE First SSR Benchmark 
System with an SVC placed at the midpoint of transmission line. The efficacy of 
various individual and composite signals is investigated at all the four critical 
levels of series compensation. The following conclusions are made : 

• Line current signal can damp all modes only at very high series 
compensation level of 90.8 %. 

• The performance of remote frequency signal is better than CIF when 
considered in combination with line current. 

• A combination of line current and remote frequency auxiliary controller 
is able to stabilize all the system modes for all the critical series 
compensation levels. Similar results were not obtained with combination 
of line current and CIF. 

• Both the remote frequency and line current signals require phase lead 
controllers. 


73 



Chapter 5 


STUDY OF TELECOM DELAYS 
AND LOCATION OF SVC ON 
SSR MITIGATION 


5.1 Introduction 

In the previous chapter, it has been studied that a combination of line 
cunent and remote frequency signal is able to damp all the torsional modes at all 
critical series compensation levels. In this chapter the influence of telecom delays 
associated with remote frequency signal on the stability of torsional modes is 
investigated. 

In this chapter an attempt is also made to study the effectiveness of 
different auxiliary signals when SVC is placed at the generator end. 

In this chapter, system studies are specifically done for series compensation 
level of 35.6 %, at which mode 4 is critically destabilized. This compensation level 
is chosen because eigenvalue results presented in Table 1.1 for IEEE First SSR 
Benchmark System reveals that the undamping of mode 4 at 35.6 % is much 
higher than undamping of modes 3 and 2 at 54.8 % and 73 % respectively. 
Although, undamping of mode 1 at 90.8 % is higher, this level of series 
compensation is too high to come in realistic range. 

5.2 Influence of Telecom Delays on System Stability 

In the previous chapter it has been studied that with a telecom delay of 
5 ms. the system is stable with (line current + RF ) auxiliary controller for all 
the critical levels of series compensation. In this section, the effect of varying 
telecom delays on system stability has been studied for critical series compensation 
level of 35.6 %, when mode 4 is critically destabilized. Eigenvalue analysis is 
done for telecom delays of 1.0 ms, 2.5 ms, 5.0 ms and 7.5 ms. Eigenvalues with 
varying telecom delays for (line current + RF) auxiliary controller are presented m 

Table 5.1 . 
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TABLE 5.1 

t td ) of remote frequency signal on system 




TABLE 5.2 

SYSTEM EIGENVALUES WHEN SVC IS PLACED AT THE GENERATOR END . 

LINE SERIES COMPENSATION = 35.6 % 
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5.3 Influence of SVC location on the effectiveness of 
different auxiliary signals 

In this section, the effect of line current and generator frequency auxiliary 
signals is examined when SVC is placed at the generator end, for critical series 
compensation level of 35 6 %, i.e. when mode 4 is critically destabilized. Table 5.2 
presents the eigenvalue results with line current and generator frequency auxiliary 
signals when SVC is placed at the generator end. Both these signals are 
considered individually. 

5.4 Discussions 

• From eigenvalue results presented in Table 5.1, it is concluded that at 
critical series compensation level of 35.6 %, the system is stable for 
telecom delay (T T d) ranging from 1.0 ms to 5.0 ms. However for 
Ttd=7.5 ms, the system was unstable. As telecom delay is reduced 
from 5.0 ms to 2.5 ms, damping of mode 4 is increased approximately 
8 times whereas that of electrical mode by 5 times. If telecom delay is 
further reduced i.e. from 5.0 ms to 1.0 ms, then damping of mode 4 is 
increased 13 times and of electrical mode by 10 times. 

Thus we see that telecom delay associated with remote frequency signal 
has significant effect on the stability of critical mode (in present case 
mode 4 ) and the electrical mode, whereas rest of the system modes are 
not much affected. The performance of (line current + remote 
frequency) auxiliary controller when tuned to modes 3 and 2, i.e. at 
critical compensation levels 54.8 % and 73 % has been investigated in 
chapter 4 with a time delay of 5.0 ms which is observed to give 
pessimistic results. If telecom delay can be reduced, then the system 
stability can be further improved. 

SVC located at Generator End 


• From eigenvalue results presented in Table 5.2, it is concluded that 
when SVC is placed at the generator end, line current auxiliary 
controller alone could not stabilize mode 4, but generator frequency 
auxiliary controller alone was able to stabilize whole system at critical 
compensation level of 35.6 %. If we compare the eigenvalue results of 
generator frequency auxiliary controller presented in Table 5.2 with that 
of (line current + remote frequency) auxiliary controller presented in 
Table 4.1, both for critical series compensation level of 35.6%, then it 
can be concluded that an SVC placed at the generator end is much 
more effective than midpoint located SVC in damping the torsional 
oscillations. This behavior is very much expected because when SVC 
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is placed at the generator end then there is no time delay associated 
with generator frequency signal, hence generator frequency signal alone 
is much more effective than composite line current + remote frequency 
auxiliary controller. 

It is important to note that if an SVC is placed at the generator end 
then it is solely dedicated to the damping of subsynchronous oscillations 
[14] and can not be simultaneously employed for any other purpose 
such as improvement of dynamic stability limit. Since SVC is an 
expensive device, therefore midpoint location of SVC is more justified 
both from technical and economical points of view. 

5.5 CONCLUSIONS 

• The telecom delay associated with remote frequency signal plays a 
significant role in system stabilization. The lower the telecom delay, 
better is the damping of torsional modes. A fast communication 
technique such as fibre optics is therefore recommended. 

• It has been revalidated that an SVC placed at generator terminals and 
employing rotor frequency auxiliary feedback is able to damp 
subsynchronous oscillations. 
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Chapter 6 
CONCLUSIONS 


6.1 General 

Static var compensators (SVC) are traditionally placed at the generator 
terminals to mitigate subsynchronous oscillations (SSO) in series compensated 
electric power transmission systems. An auxiliary control [14] of the SVC 
employing generator rotor frequency feedback is shown to be adequate for 
damping SSO. These SVCs suffer from the disadvantage that they are exclusively 
employed for SSO suppression and can not be utilized for achieving any other 
objective simultaneously. It may be understood that SVCs are very expensive 
equipment which need to be utilized most optimally. 

It has been reported [15] that an SVC located at the midpoint of a series 
compensated line can be utilized for dual purposes of damping SSO and stability 
enhancement. Controllers of these SVCs have been designed based on a 
combination of line current signal and computed value of generator internal 
frequency (CIF). Both these signals utilize local measurements, as the current 
trend in power systems practice is to use local signals for the sake of reliability. 

In this thesis an altogether new concept of SVC control employing a 
remote generator frequency signal transmitted over a telecom line, is presented. 
The motivation for this new investigation is provided because of two reasons : 

(i) Signals such as “current margin” are successfully transmitted over 
telecom lines in HVDC transmission systems in a reliable manner, 
over the world. 

(ii) Since generator rotor frequency is the most effective signal when 
SVC is connected at the generator terminals, it is thought that the 
same signal, if made available at the remotely located SVC, may 
also prove to be substantially effective. While the previous work 
[15,23] attempted to synthesize this signal based on the local 
measurements, in this thesis the same rotor frequency signal is 
actually measured and transmitted over telecom channels. 


79 



In this thesis the study system is chosen to be the IEEE First SSR 
Benchmark System. This is suitably modified to connect an SVC at the midpoint 
of the transmission line. Different auxiliary signals such as line current, CIF, 
remote frequency signal and combination of line current with CIF and remote 
frequency signal are examined in respect of their effectiveness in stabilizing the 
torsional oscillation modes. The auxiliary controller is, a priori, assumed to be a 
simple lead - lag transfer function. System studies are done for all the four critical 
series compensation levels, i.e. 35.6%, 54.8%, 73% and 90.8% at which mode 4, 
mode 3, mode 2 and model, respectively, are critically destabilized. 

It is observed that a pure voltage control of SVC has minimal effect on 
the stability of torsional oscillation modes, which is expected since voltage control 
does not contribute electric damping. When the network is critically compensated, 
which may not be realistic, a combination of line current and remote frequency 
signal is successful in damping torsional oscillations at critical compensation levels 
of 35.6 %, 54.8 % and 73 %. At critical compensation level of 90.8 %, line 
current auxiliary controller alone is able to stabilize the system. Eigenvalue results 
of this thesis clearly show that (line current + RF) combination is much superior to 
(line current + CIF) even with a time delay of 5 ms. 

It has been shown that with a time delay of 5 ms, the study system is 
stable with (line current + remote frequency ) auxiliary controller over whole range 
of critical series compensation. It is further shown that the lower is the value of 
telecom delay, better is the stability of torsional modes. 

It is concluded in this thesis that this remote rotor frequency signal 
together with line current signal can be successfully applied for damping all the 
torsional modes at all the critical levels of series compensation. The performance 
of remote rotor frequency signal is much better than the computed internal 
frequency signal. 

This concept of remotely transmitted signal is being investigated for the 
first time in the control of static var compensators to damp subsynchronous 
oscillations. 

6.2 Scope for Future Work 

• To study the capability of SVC auxiliary controllers in damping 
torsional oscillations for IEEE Second Benchmark System for SSR 
and even on big sized realistic systems . 

• To study the performance with exact CIF auxiliary controller for 
midpoint located SVC , where transients associated with inductor and 
capacitor are taken into account . 
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• To validate the eigenvalue results obtained in chapter 4 and 5 through 
EMTDC/PSCAD simulation , in which it is possible to model all the 
system nonlinearities . 
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APPENDIX A 


SYNCHRONOUS MACHINE MODEL 

PARAMETERS 

In this Appendix, expressions are given for the various constants which are 
used in the synchronous machine model [6] . 

The constants a, - a 9 are defined as 


'a, a 2 

W o 

1 

'*'s 

* 

> 

Rj x j}, 

_ a i «4. 

XjX.-X 2 ^ 

R h X Jh 

R h Xf 


A 

<*6 

W o 

’ R g X k 

' ~ R S V 

.°1 

a s_ 

X g X k ~ * A 

_~ R k X * 

** 

1 


a 9 = -1 / T CO(/ ; T C0I , is the time constant of dummy coil 

(small value chosen arbitrarily = 0. 1 msec ) 

The constants b x - b 6 are defined as 


' V 0 R f 

i 

O- 

1 

W o 

R /( X df X h X dh X fh ) 

X df 

M 

X f X„-X 2 fl, 

R h( X / X dh ~~ X Jh X (lf )_ 


V 

W o 

X< x « x * - Vg *>" 


A. 

x g x k -X 2 gk 

R k( x g x <,x - 


^4 — ^ X c^ X q X d ) 

w 0 = 2nf ; / is the system frequency 
The constants c, - c 5 are given as 
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X d/ X h ~ x dh x Jh 
X d\ x f x h- x Jh 2 ) 


_ X dh X f X jh X df 

x 7 ^r^ 7 ~) 


_ y k-Vgk 

x d\ x g x k- x gk 2 ) 

__ ViZVsi 

x d\ x g x k~ x gk 2 ) 



where 

x f ,x h ,x g ,x k are reactances of the rotor coils specified by the 
subscripts. 

x c is the reactance of the dummy coil (chosen = x q ) 

Ry , R h , R g , R k are resistances of the rotor coils specified by the 
subscripts. 

x df , x dh , Xjb , x gk , x qg , x qk are mutual reactances between rotor coils 
specified by the subscripts. 

The resistances and reactances of the various rotor coils are defined as 
follows : 

X q g = X qk = X qc = X gk = x q ~ x i 

(x d ~ x i) x df 
X/j - 7 

( x d ~ x d ) 


x df ~ X dh ~ X jh X d X l ’ 


X 


( x d~ X d) 




(V -*/)(*, -X,) 

( x o ~ X q) 


X kl = 
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X f =X df+ X JI 


X h X df "*■ X hl 


x g = x qk+Xgl 


X k = X qk + X kl 


R, 


l df 


W o T dA x d ~ x d ) 


R, 


_ zi£. 

W o T do\ x d ~ x d) 


R n 


(Xq -X ,) 1 


W o T qo ( X q ~ X q ) 




V 


where 


W o T qo ( X q~ X q ) 


X/ is the stator leakage reactance. 

x d , x d , x d are the direct axis synchronous , transient and subtransient 
reactances respectively. 

x q , x q , x q are the quadrature axis synchronous , transient and subransient 
reactances respectively. 

T do > T do ' are the direct axis transient and subtransient open circuit time 
constants respectively. 

T q0 , T qo are the quadrature axis transient and subtransient open circuit 
time constants respectively. 
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APPENDIX B 


DETAILS OF SYSTEM MODEL 
DESCRIBED IN CHAPTER 2 


B.l Rotor Circuits 


Substituting eqn.(2.7) in eqn.(2.6) 

4^ = a^ f +a 2 ^ h + b x v f + b 2 cos 8 i D - b 2 sin£ i Q 
y ¥ h =a^ f + a 4 + b 3 cos S i n ~ b 3 sin §i Q 

'Pjj = a 5 x i J g +a 6 y ¥ k +b 5 sin Si n +b 5 cos Si Q (B.l) 

x ¥ k = a 1 x ¥ g + a z '¥ k + b 6 sin Si D + b 6 cos Si Q 
l P c = a g y f c + b 4 sin Si D + b 4 cos Si Q 


The state equation for the rotor circuits is obtained by linearizing eqn. (B.l) 


Xr ~~ ArXr + Br\Ur\+BriUri 
wherej- 

x, = K A4/ /. A ^ A ^ A ^ } . UrA^ M 

UR2 = l Ai » A/ J' 



a 2 

0 

0 

0 " 


~^2 

O' 

a 3 

a 4 

0 

0 

0 


~h iqO 

0 

0 

0 

a s 

a 6 

0 

’ Br\ 

V* o 

0 

0 

0 

«7 

a s 

0 


^6 idQ 

0 

_0 

0 

0 

0 

a 9_ 


. h'dt) 

0 


(B-2) 
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B 


R2 


b 2 cosA 0 

~b 2 sinS 0 

b 2 cosS 0 

- b 3 sin A 0 

b s sinA 0 

b } cosS 0 

b(, sin Sq 

b 6 cos A 0 

b A sin A 0 

b 4 cos A 0 


ido ~ cos ^o‘do ~ o^qo 

l qO = sin <5 (ji qq + cos S 0 i g0 

Substituting eqn. (2.9) in eqn. (2. 1 0) gives 

I D =c x cosS y V f + c 2 cosA v F / , + c 3 sinA4' g + c 4 sinA'f' t + c 5 sinA l F c 

(B-3 ) 

I Q = -c, sinATy - c 2 sin A x ¥ h + c 3 cosAT^ + c 4 cosA v J( t + c 5 cosA v F c 


Differentiating eqn. (2.9) 


/ < ,=c 1 4' / +c a V* 

I^C 3 %+C 4 %+C 5 % 

Eqn. (2.10) is also differentiated to give 


cos A 

sin A 

"id 

dS 
+ — 

-sin A 

cos A 

-sin A 

cos A 

K. 

dt 

-cos A 

-sin A 


...(B.4) 


(B.5) 


The flux derivative terms in eqn. (B.4) are replaced by their respective 
expressions obtained from the rotor state equations (B.2 ) . I d , I q so obtained 
are substituted in eqn . (B.5) . The resulting equation together with eqn. (B.3) are 
now linearized to give the output equations of rotor circuit as 


y R rCmX R + DmUm (B-6) 

y R2 ~C^X R + D R 2U K ^DmU R 2 (B-7) 

where = [a I d A I Q ] 1 ,y ^ = [a 7 d A I Q ] ' 

The nonzero elements of the different matrices are defined as 


C RX (1,1) = c, cosA 0 
C R\ 0 >2) = c 2 cos Sq 
C R i (1 ,3) = c 3 sin A 0 
C R] (1 ,4) = c 4 sin A 0 


86 



C RX (1,5) = c 5 sin S 0 
C R[ (2,1) = -c, sinS 0 
C/?i (2,2) = -c 2 sin£ 0 
Cr\ (2,3) = c 3 cos 
C R1 (2 ,4) = c 4 cos £ 0 
C*. (2,5) = Cj cos^ 0 

C R2 (l,l) = (^l«l + C 2 flr 3) COS ^0 
C R2 (1,2) = (c x a 2 + c 2 a 4 ) cos 5 0 
C R2 (1 3) = (c 3 a 5 + c 4 a 7 )sin S 0 
C R2 (1,4) = (c 3 a 6 + c 4 a 8 )sin <5 0 
C R2 (1,5) = c 5 a 9 sin^ 0 
C R2 (.2,\) = -(c x a x + c 2 a 3 )sinS 0 
Q 2 (2 ,2) = -(c,a 2 + c 2 a 4 ) sin 5 0 
C R2 (2,3) = (c 3 a 5 +c 4 a 7 )cos(5 0 
Crz (2,4) = (c 3 a 6 + c 4 a 8 )cos<5 0 
C R2 (2,5) = c 5 a 9 cos 5 0 
Dri (1 4) = Iqo 
D ri (2 ,1) = - I do 

D r 2 (1 ,\) = -(c l b 2 + c 2 b 3 )i q 0 cosS 0 + (c 3 b 5 + c 4 b 6 + c 5 b 4 )i d 0 smS 0 


D R2 {l,2)=I Q o 

D r 2 (2 ,1) = (c,& 2 + C 2^3 ) ' 9 0 sin ^0 + ( C i b 5 + C 4 & 6 + C 5^4 ) J rfO COS S 0 


D r 2 (2 ,2) = - / D0 
£> fi3 (1,1) = c,&,cos£ 0 

D R3 (2,l) = -c I 6 1 sin^ 0 

D ft4 (1 ,l) = (C,f>2+C 2 & 3 )COS 2 <^0 + (c 3 ^5 +C 4 6 6 +C 5 6 4) sin ^0 

D R4 (\ ,2) = ~{c { b 2 + c 2 b 2 ) sin 5 O cos S 0 + (c 3 b 5 + c 4 b 6 + c 5 b 4 ) sm S 0 cos <5 0 

D R4 (2 ,1) = -(c,& 2 + c 2 6 3 ) sin 5 O cos 6 0 + {c 2 b 5 ■ + c 4 b 6 + c 5 b 4 ) sin <5 0 cos 8 0 
D r4 ( 2 , 2 ) = (c x b 2 +c 2 b 2 )sm 2 S 0 + (c 3 b 5 + c 4 b 6 + c 5 b 4 ) cos <5 0 


I D0 =cos<5 0 / d0 +sin<V 9 <) 

7g 0 = -sin^g/jo +cos<5 0 / 9 q 


B.2 State Space Model of the Turbine - Generator Mechanical 
System 

The six mass representation of the turbogenerator shaft system is modeled 
in Fig B 1 wh”“ M represents the moment of inertia , K the spnng constant , 
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D is the damping coefficient and 8 represents the angular position of 
concerned mass . T mi is the mechanical torque input corresponding to i th mass 
while T e is the electrical torque on the generator. Equations relating the angular 
position and velocity of various masses are given as 

+ (£>„ + D u )S l -D n S 2 + K ll S l -K l2 S 2 = T ml 

M 2 S 2 ~ E)|2<£ |+ (£> )2 + D 22 + D 22 )(5 2 — D 22 8, — K,jS, 

+ (K u +K 2i )8 2 -K 2} 8 3 =T m2 

^38,— D 23 S 2 + (D 2} + £) 33 + D u )S 3 — Dj 4 S 4 - K 23 5 2 

+ (K 23+ K 3i ) 5 3 -K 34 5 A =T mi 

MA — D 3A 5 3 + (£>34 + D AA + D\s )S 4 — D as S 5 — K 34 S 3 

+ (K 34 +K 45 )S A -K 45 S 5 =T mi 

M 5 S 5 - Dj A + (£> 45 + £> 55 + £> 56 )S 5 - D 56 4 - K A5 S A 

+ {K 45 +K 56 )8 5 -K 56 S 6 =-T e 

M 6 S 6 -Dj 5 + (£> 56 + £> 66 )S 6 - K 56 8 5 + K 56 S 6 = 0 
Linearizing eqns. (B.8) and writing in state space form as 


Xu Am Xm ~ >r Bm Um 


(B.9) 


where 

Xu =[a 5 , a s 2 as, as a as 5 as 6 as, 

Um = [A T ml A T m2 A £ m3 A £ m4 -A T e ] 


A S 2 AS 3 a 4 a 4 a 4J 


Am 


0 I 
Ami Ami. 


I is an identity matrix of dimension ( 6 , 6 ) . 


The nonzero elements of Ami ^ defined by 


A„(U)=- 


*!2 

Mj 


Ami(^ = 


£2 
M ( x 
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K, 


A( 2,1) = ^ 

AmK A/, 


(^12+^23) 


Au\( 2 ^ 


K 


23 


m 7 


K 


Am( 3 ’ 2 y 

(K 23 + -^34) 


Arn^y 


M-, 


Am (3>4) 




34 

Mt 


. (^34 +*45) 

^> (4,4) m 4 
^w' (4,5)= m7 

A r S) _ (*45 +*56). 

u4a/i C 5 ^ 5 ) 


^><l (6 ’ 5) M 


6 

*56 


^a/i^ 6 ’ 6) M 6 
The nonzero elements of y4 w2 are § iven 


as 
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Am 2 AO 
AmiM 
Am 2 ( 2,1) 
AwS 2 ’ 2 ) 
Am < 2 ’ 3 ) 

Am 2 C> 2) 
^, 2 ( 3 >3): 
v4„ 2 ( 3 > 4 ) 

a, 2 ( 4 >3)= 

^ W2 ( 4 ’ 4 ) 

A, 2 ( 4 ’ 5 ) 

Am 2 ^ 

AmW 

Am 2 ^’ 6 ) 


(A.+A?) 

K 


A 2 

M x 

£n 
M 2 


(Ai + A 2 + A3) 
m 2 


As 

m 2 

As 

m 3 


(P 23 + £) 33 4 - £>34) 

m 3 


Ai 

m 3 

Ai 

m 4 


(Aj + Ai + As) 
a/ 4 




(As + As + As) 
a / 5 


Ak 

m 5 


91 





M 6 


,<W 6 > 6 ) 


A 


66 


M, 


Since self damping and mutual damping between the masses is assumed to 
be zero , for IEEE First Benchmark System . Matrix Am 2 becomes a zero 
matrix . 


Am 


0 / 

Am 0 


Since the governor system involves large time constants , it is ignored in 
the present analysis . Thus , 


A T m{ = A T ml = A r m3 = A7 m4 = 0 


(B.10) 


The electromagnetic torque T e acting on the generator rotor, as given 
in [6 ] , is expressed by 


T e = -x d ( i d I q —i q I d ) (B. 11) 

The currents i d ,i q and I d I q are transformed to D-Q reference frame 
using the relationship given in eqn. (2.7) . It is noted that , though eqn . (2.7) is 
written specifically for transforming i d , i q to D -Q axis quantities , the same 

relation also applies for the transformation of I d I q to corresponding currents in 
D-Q frame of reference as 


T e =-x d [(cos<? i D - sin 5 iQ)(smS I D +cos 8 1 q) 

- (sin 8 i D + cos 8 i Q )(cos 8 I D - sin 8 Iq )] 

Alternatively , T e = -x d ( i D lg - ‘q^d) 
linearizing above eqn . we get , 

A T e = —x d [A + 'do A A ~ A iglpo ~ *goA-A>] 

Substituting eqns. (B.10-B.14) in eqn . (B.9) results in 

X M = Am Xm + Bm\ U m + Bm2 Umi 
where 


(B.12) 

(B.13) 

(B.14) 


(B.15) 
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B m { 11,D= 

B An(H,2) = 

5a/ 2 0U) = 


A/p]' 

U M2 

[ a *d 

A ip 

t 

elements of matrices 

Bm 1 

and 

B 

x d ’qo 





M 5 





: d ’DO 






5 a/2 (11,2)=- 


M 5 


X d IqO 


M 5 

A/5 


The output equation of the mechanical system is given by 

(B.16) 

y t rCM 

where 

y u =\ AS * AA >1 

The nonzero elements of Cm ^ given by 
C a/ (1,5) = 1-0 

C w (2,ll) = l-0 
It may be noted that 

AS = AS 5 (B.17) 

Aco = A <5 5 

where A<5 and Ant are the incremental changes in generator rotor 
angle and angular speed respectively . 

B.3 Network 

The nonzero elements of the matrices utilized in eqn. (2.19) are given 

below 
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5,(U) = 


R, 


S\ (1»2) = — 7- 


51(24) = 


1 


5,(14) 


= _ 


5 3 (M) = 


The different matrices in the state equation of the network given by 
eqn (2.22) are defined as follows 


A, 


O) 0 I 


*oi Si 


B, t = 


0 ® 0S2 

-®. S 2 0 


B A 2 ~ 


S 2 0 

.° S 2 _ 

The output equation for the network model is written as 


y N C N X N 

where = [a/ d M q ] ' 

The nonzero elements of Cn are defined by 
C N (1,1) = C N (2,3) = 1.0 

B.4 Interconnection of Various Subsystems 

The various vectors and matrices are defined as follows 

Xz = X R Xaz xj' 

Ur\ Ur 2 U\f\ Um2 Un\ Uni\ 

) 


Ui = . 

y,=\y 


A 


(B.i8) 


ri y R2 y m y n 
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Bm Bn 0 0 0 0 

Br= 0 0 B w B\n 0 0 

0 0 0 0 #v, 


C* 0 0 

r _ C« 0 0 

0 C, 0 

o o c N 


Dr 


Dm 

0 

0 

0 

0 

0 

Dm 

Dm 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 
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APPENDIX C 


DETAILS OF SYSTEM MODEL 
DESCRIBED IN CHAPTER 3 


C.l Network Model 

The nonzero elements of matrices Si , S 2 , S 3 and S 4 described for 
the network in eqns . (3.6) , (3.7) are defined as 


5 , 0 , 1 ) = -^ 

5, (1,4) = - 7 - 

l b 

S\ (2,3) = -j— 
S,(3.1) = ~ 
S,(4. l) = d- 

^se2 

s 2 a D=-d- 
s,ax)=-~- 

l a 


$i (1>3) = — ■ 

S i (2,2) = ~~ l 

^A 

S,(2,5) = ~ 

L A 

5,(3,2) = 
S,(5,2) = -J- 


s.(U)=-d- 

The different matrices in the state equation of the network given by eqn . 
( 3 . 8 ) are defined as follows : 



The 


„ . s, --• r 

As \o)j Si 


Bsi 


0 

-6> a Si ° 


rs, 0 

Bs-[ o S, 

fS» 0 

Bs 3= 0 s 


(C.l) 

(C.2) 


. „ f the network model are given as 
tput equations of 

y M =CxiX.v 

V„ S Cn 2^ r V 

i y =[AV 3 D AV 3q1 

where A '° ’ ' d Cw are defined by 

r (U) = C«. (2,7)= ' ° 

* static VAR System ^ used * C^> “ 

T* nonzero elements of the 


X(U) = 


CO o 


Q 

A,(2,l) = «o 

T s C3,4) = ^-° 

K p fipo 

K, 

A. s (5,3) = ~ r 


i s (2,2) 


cy o 
<2 


is( l,6) = n>o v 3»o 
j 4 s ( 2,6) = ® o v 3a° 


As- (6,5) : 


i s (4,2) = 
/ 4 s (5,4) = 

A s (6,6) 


K i) hg o 
Tf.t ho 


Kp 

T s 


,As(4>4) = 

A s ( 5 , 5 ) : 


l 

'T, 


M 


‘Ts 


r D 


>0 
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S 5 ,(l,l) = ffl„S 0 


B sl ( 2 , 2 ) — (o 0 Bq 


B, ( 4 , 1 ) = -^ . B S] ( 4 , 2 ) = —~— 

V 30* M v 

S„(3,l) = 1.0 , S s! (5,1) = -|l 

*s 

C. si (14) = Cy, (2,2) = 1 .0 , £) s = 0 

C.3 Interconnection of Various Subsystems 


In this case , the various vectors and matrices are defined as follows : 
X 7 = X R Xu Xn Xsi 

Hi ~ 11 r\ Ur 2 Mu i M\n Ms\ M.\z. Msi W.vJ 

y,=\y„ y„ y„ y>, y«, y) 

'0 0 / 0 0 0 
0 0 0 / 0 0 

/ 0 0 0 0 0 

0 0 0 / 0 0 

*> 0 0 0 0 0 / 

/ 0 0 0 0 0 

0 / 0 0 0 0 

0 0 0 0 /0 

~Bn B R 2 0 0 0 0 0 °~ 

D = 0 0 Bux Bm 0 0 0 0 

0 0 0 0 Bn Bsi Bm 0 

0 0 0 0 0 0 0 Bs 1. 

c o 0 o’ 

\-sR\ 

r o o o 

\~sR2 

r 0 c„ 0 0 

Lr 0 0 c»-, 0 

o 0 Cm 0 

o 0 0 Cs. 


A, 0 0 0 

. 0 A,, 0 0 

A o » A » 

ooo As 
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D t = 


Dm 

Dm 

0 

0 

0 

0 


0 0 0 0 0 0 0 
Dn 000000 
0 0 0 0 0 0 0 
0 0 0 0 0 0 0 
0 0 0 0 0 0 0 
0 0 0 0 0 0 0 


C.4 Details of Auxiliary Controller Model 

C.4.1 Line Current Auxiliary Controller 

Eqn. (3.32) can be written as 

"r=F,Xr (C.3) 

where u ( = A i 4 

The nonzero elements of p are given by 
^(1,19) = --^ , TJ(1,24) = — -^1 

^40 *40 

Eqn. (C.3) can be rewritten in matrix form as 

u c — P (R Xr + FcuXu + FcnXn + FcsXs (C-4) 


where the various matrices are defined as submatrices of p } 

If„ Fcu F„ Fj=F, 

It is noted that f ( R = f cu = Fes = 0 
The nonzero elements of p cw are given by 

F w (1,2) = F i (U9) 

F cn { 1,7) = F 1 (1,24) 

The state and output equations of the general auxiliary controller depicted 
in Fig. 3.8 are written as 

Xc = AcX< + B( u < (C 

y<=CcXc + Dc u c (C 


99 



where X( = [z 5 ] 


Jr = Av, 


The nonzero element in each of the matrices y£ c , J$ ( . , (J c and £) 
are defined as 

A(U)=~ 

1 2 

h **2 

C, (1,1) = 1.0 

A(U) = A:»i 

h 

C.4.2 Computed Internal Frequency (CIF) Auxiliary Controller 

Linearizing eqns. (3.41) , (3.42) and substituting in eqn. (3.44) gives 

M,= F,Xr + F,Xr < C7 > 

where the nonzero elements of p' 3 and J7 A are defined as 

F 3 (l,l9) = co 0 L E ^ , F 3 (l,20) = -% 

e 0 e o 

F 3 (l,24) = o 0 L E ^ , F 3 ( 1,25) = % 

e 0 e 0 

F 4 (U9) = T,% , F, (1,24) = -L F 

e o e o 

As the perturbation in SVC reference voltage , AV k/ , is not considered in 
the analysis of auxiliary SVC control as presented in chapter 3 , eqn.(3.30) gets 
modified to 

Xt=AXt + B u s3 (C8) 

where u S3 = Av ; 

Matrices y£ and are already defined in chapter 3 . 

Substituting eqn. (C.8) in eqn. (C.7) gives 

AS i = F 3 Xr + F SA Xi + B u : sj) 

or 
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(C.9) 


MrFsXr 

where 

F-= F, + F,A • F,B=o 

The CIF signal is obtained by differentiating eqn. (C.9) 
Aco i = F 5 X 7 

Resubstituting % T from eqn.(C.8) in eqn. (C. 10) results in 
A g>i=F*(A Xt + B u si ) 

a ®^F 6 Xt 

where F 6 = FsA > 

The nature of matrices F 5 and ft ls suc ^ ^ at 

F,B = <> 

Eqn. (C. 11) can now be written in expanded form as 

u c = F cr Xr + F cm Xm + F cnXn + FcsXs 
where u c = Aco , 

The various matrices are defined as submatrices of 33 

L Fcr Fcm Fcn Fa\ = Fi 


(C.10) 


(C.ll) 


(C.12) 
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APPENDIX D 


SYSTEM DATA FOR IEEE FIRST 
SSR BENCHMARK MODEL SYSTEM 

D.l System Base Quantities 

Base Voltage = 500 kV 
Base MV A = 892.4 
Base Frequency = 60 Hz 

D.2 Generator Data 

Sn = 892.4 MV A 
Ra = 0.0 pu 

x d =\.19 pu 
x ? =1.71 pu 
x,=0.13 pu 

x' d = 0.169 pu 
x; =0.228 pu 

x^ =0.1 35 pu 
x; =0.200 pu 

T!, o=4.3 sec. 

Ko= 0.85 sec. 

r; 0 = 0.032 sec. 
r; o =0.05 sec. 
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D.3 Multimass Parameters 


Mass 

Shaft 

Inertia 

Spring Constant 

Self 



H (sec ) 

K (pu Torque/ 

Damping 




rad ) 

D (pu) 

HP 

HP -IP 

0.092897 

19.303 

0.0 

IP 

IP - LPA 

0.155589 

34.929 

0.0 

LPA 

LPA - LPB 

0.858670 

52.038 

0.0 

LPB 

LPB - GEN 

0.884215 

70.858 

0.0 

GEN 

GEN - EXC 

0.868495 

2.822 

0.0 

EXC 


0.0342165 


0.0 


The mutual damping between the masses is assumed to be zero. The 
steady state mechanical torque is apportioned among the turbine sections HP, IP, 
LPA and LPB respectively as follows : 30 %, 26 %, 22 % and 22 %. The exciter 
steady state torque is assumed to be zero. 

D.4 Transformer Data (on system base) 

X T =Q.Upu 

D.5 Transmission Line Data 

Resistance ( R) = 0.02 pu 
Inductve Reactance {X t ) = 0.5 pu 
Capacitive Reactance (X c ) = 0.371 pu 
F ault Reactance X F ,(L-G) = 0.04 pu 

D.6 Static Var Compensator 

T m = 2.4 ms , T s =5.0 ms , Tor 5 1.667 ms 
Ki = 90 , Kp = 2.6 , K d = 0.03 pu (1%) 
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APPENDIX E 


CALCULATION OF INITIAL 
CONDITIONS 


E.l Generator Initial Conditions 

The vector diagram of the overall system is depicted in Fig. E.l .The 
generator current I g is given by 


7 


K 


* 

V , r 


(E.1) 


where V g is the conjugate of generator terminal voltage V g . 

The field voltage V f and rotor angle 8 are computed as follows 


E<", = Y* +1 A R a+J x q ) 

8 = ZE qa -n/2 


(E.2) 

(E.3) 


V,=\E qa \ +I gd (x d -x q ) 


(E.4) 


where 7^ is the component of generator current along d axis 
The initial values of other variables are 



(E.5) 
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X d^gd + X JhI f 

^ = X qg l S1 (E.6) 

8 = 0 

where I gq is the q axis component of the generator current. 
x / , x d/ -> x fh > x gq and x qk are defined in Appendix A. 


E.2 SVC Initial Conditions 

The reactive power generated by SVC fixed capacitor (FC) is given by 

Qc ~ 03 oC /( v 30 2 (E.7) 

where v 30 = initial magnitude of SVS bus voltage v 3 . 

The initial value of TCR susceptance is then computed as 

S„“(2r-fi)W) (E.8) 


where Q 3 is the net reactive power at SVC bus obtained from load flow . 
The initial value of TCR current is evaluated as 


*30 — B 0 V 3Q 


(E.9) 
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